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Abstract 

The  main  goal  of  this  paper  is  to  prove  that  if  the  energy-momentum 
(or  energy-Casimir)  method  predicts  formal  instability  of  a  relative  equi¬ 
librium  in  a  Hamiltonian  system  with  symmetry,  then  with  the  addition 
of  dissipation,  the  relative  equilibrium  becomes  spectrally  and  hence  lin¬ 
early  and  nonlinearly  unstable.  The  energy-momentum  method  assumes 
that  one  is  in  the  context  of  a  mechanical  system  with  a  given  symmetry 
group.  Our  result  assumes  that  the  dissipation  chosen  does  not  destroy 
the  conservation  law  associated  with  the  given  symmetry  group — thus,  we 
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consider  internal  dissipation.  This  also  includes  the  special  case  of  sys¬ 
tems  with  no  symmetry  and  ordinary  equilibria.  The  theorem  is  proved 
by  combining  the  techniques  of  Chetaev,  who  proved  instability  theorems 
using  a  special  Chetaev-Lyapunov  function,  with  those  of  Hahn,  which 
enable  one  to  strengthen  the  Chetaev  results  from  Lyapunov  instability 
to  spectral  instability.  The  main  achievement  is  to  strengthen  Chetaev’s 
methods  to  the  context  of  the  block  diagonalization  version  of  the  energy 
momentum  method  given  by  Lewis,  Marsden,  Posbergh,  and  Simo.  How¬ 
ever,  we  also  give  the  eigenvalue  movement  formulae  of  Krein,  MacKay 
and  others  both  in  general  and  adapted  to  the  context  of  the  normal  form 
of  the  linearized  equations  given  by  the  block  diagoanl  form,  as  provided 
by  the  energy-momentum  method.  A  number  of  specific  examples,  such  as 
the  rigid  body  with  internal  rotors,  are  provided  to  illustrate  the  results. 


1  Introduction 

A  central  and  time  honored  problem  in  mechanics  is  the  determination  of  the 
stability  of  equilibria  and  relative  equilibria  of  Hamiltonian  systems.  Of  par¬ 
ticular  interest  are  the  relative  equilibria  of  simple  mechanical  systems  with 
symmetry,  that  is,  Lagrangian  or  Hamiltonian  systems  with  energy  of  the  form 
kinetic  plus  potential  energy,  and  that  are  invariant  under  the  canonical  action 
of  a  group.  Relative  equilibria  of  such  systems  are  solutions  whose  dynamic  orbit 
coincides  with  a  one  parameter  group  orbit.  When  there  is  no  group  present,  we 
have  an  equilibrium  in  the  usual  sense  with  zero  velocity;  a  relative  equilibrium, 
however  can  have  nonzero  velocity.  Wrhen  the  group  is  the  rotation  group,  a 
relative  equilibrium  is  a  uniformly  rotating  state. 

The  analysis  of  the  stability  of  relative  equilibria  has  a  distinguished  history 
and  includes  the  stability  of  a  rigid  body  rotating  about  one  of  its  principal  axes, 
the  stability  of  rotating  gravitating  fluid  masses  and  other  rotating  systems. 
(See  for  example,  Riemann  [1860],  Routh  [1877],  Poincare  [1885,  1892,  1901], 
and  Chandrasekhar  [1977]). 

Recently,  two  distinct  but  related  systematic  methods  have  been  developed 
to  analyze  the  stability  of  the  relative  equilibria  of  Hamiltonian  systems.  The 
first,  the  energy- Casimir  method,  goes  back  to  Arnold  [1966]  and  was  developed 
and  formalized  in  Holm,  Marsden,  Ratiu,  and  Weinstein  [1985],  Krishnaprasad 
and  Marsden  [1987]  and  related  papers.  While  the  analysis  in  this  method  often 
takes  place  in  a  linear  Poisson  reduced  space,  often  the  “body  frame” ,  and  this  is 
sometimes  convenient,  the  method  has  a  serious  defect  in  that  a  lack  of  sufficient 
Casimir  functions  makes  it  inapplicable  to  examples  such  as  geometrically  exact 
rods,  three  dimensional  ideal  fluid  mechanics,  and  some  plasma  systems. 

This  deficiency  was  overcome  in  a  series  of  papers  developing  and  applying 
the  energy-momentum  method;  see  Marsden,  Simo,  Lewis  and  Posbergh  [1989], 
Simo,  Posbergh  and  Marsden  [1990,  1991],  Lewis  and  Simo  [1990],  Simo,  Lewis, 
and  Marsden  [1991],  Lewis  [1992],  and  Wang  and  Krishnaprasad  [1992].  These 
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techniques  are  based  on  the  use  of  the  Hamiltonian  plus  a  conserved  quantity. 
In  the  energy-momentum  method,  the  relevant  combination  is  the  augmented 
Hamiltonian.  One  can  think  of  the  energy  momentum  method  as  a  synthesis  of 
the  ideas  of  Arnold  for  the  group  variables,  and  those  of  Routh  and  Smale  for  the 
internal  variables.  In  fact,  one  of  the  bonuses  of  the  method  is  the  appearance 
of  normal  forms  for  the  energy  and  the  symplectic  structure,  which  makes  the 
method  particularly  powerful  in  applications. 

The  above  techniques  are  designed  for  conservative  systems.  For  these  sys¬ 
tems,  but  especially  for  dissipative  systems,  spectral  methods  pioneered  by  Lya¬ 
punov  have  also  been  powerful.  In  what  follows,  we  shall  elaborate  on  the 
relation  with  the  above  energy  methods. 

The  key  question  that  we  address  in  this  paper  is:  if  the  energy  momentum 
method  predicts  formal  instability,  i.e.,  if  the  augmented  energy  has  a  critical 
point  at  which  the  second  variation  is  not  positive  definite,  is  the  system  in  some 
sense  unstable ?  Such  a  result  would  demonstrate  that  the  energy-momentum 
method  gives  sharp  results.  The  main  result  of  this  paper  is  that  this  is  indeed 
true  if  small  dissipation,  arising  from  a  Rayleigh  dissipation  function,  is  added 
to  the  internal  variables  of  a  system.  (Dissipation  in  the  rotational  variables 
will  be  considered  in  another  publication.)  In  other  words,  we  prove  that 

If  a  relative  equilibrium  of  a  Hamiltonian  system  with  symmetry  is 
formally  unstable  by  the  energy-momentum  method,  then  it  is  lin¬ 
early  and  nonlinearly  unstable  when  a  small  amount  of  damping 
(dissipation)  is  added  to  the  system. 

Some  special  cases  of,  and  commentaries  about,  the  topics  of  the  present 
paper  were  previously  known.  As  we  shall  discuss  below,  one  of  the  main  early 
references  for  this  topic  is  Chetaev  [1961]  and  some  results  were  already  known 
to  Thomson  and  Tait  [1912].  See  also  Ziegler  [1956],  Haller  [1992],  and  Sri 
Namachchivaya  and  Ariaratnam  [1985].  The  latter  paper  shows  the  effect  of 
dissipation  induced  instabilities  for  rotating  shafts,  and  contains  a  number  of 
other  interesting  references. 

Next,  we  outline  how  the  program  of  the  present  paper  is  carried  out.  To 
do  so,  we  first  look  at  the  case  of  ordinary  equilibria.  Specifically,  consider  an 
equilibrium  point  ze  of  a  Hamiltonian  vector  field  Xh  on  a  symplectic  manifold 
P,  so  that  Xj{(ze )  =  0  and  H  has  a  critical  point  at  ze.  Then  the  two  standard 
methodologies  for  studying  stability  mentioned  above  are  as  follows: 

(a)  Energetics  —  determine  if 

62H(Z')  =:  Q 

is  a  definite  quadratic  form  (the  Lagrange- Dirichlet  criterion). 
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(b)  Spectral  methods  —  determine  if  the  spectrum  of  the  linearized  op¬ 
erator 

D XH(zt)  =:  L 

is  on  the  imaginary  axis. 

The  energetics  method  can,  via  ideas  from  reduction,  be  applied  to  relative 
equilibria  too  and  this  is  the  basis  of  the  energy-momentum  method  alluded  to 
above  and  which  we  shall  detail  in  §3.  The  spectral  method  can  also  be  applied 
to  relative  equilibria  since  under  reduction,  a  relative  equilibrium  becomes  an 
equilibrium. 

For  general  (not  necessarily  Hamiltonian)  vector  fields,  the  classical  Lya¬ 
punov  theorem  states  that  if  the  spectrum  of  the  linearized  equations  lies  strictly 
in  the  left  half  plane,  then  the  equilibrium  is  stable  and  even  asymptotically  sta¬ 
ble  (trajectories  starting  close  to  the  equilibrium  converge  to  it  exponentially  as 
t  — ►  oo).  Also,  if  any  eigenvalue  is  in  the  strict  right  half  plane,  the  equilibrium 
is  unstable.  This  result,  however,  cannot  apply  to  the  purely  Hamiltonian  case 
since  the  spectrum  of  L  is  invariant  under  reflection  in  the  real  and  imaginary 
coordinate  axes.  Thus,  the  only  possible  spectral  configuration  for  a  stable  point 
of  a  Hamiltonian  system  is  if  the  spectrum  is  on  the  imaginary  axis. 

The  relation  between  (a)  and  (b)  is,  in  general,  complicated,  but  one  can 
make  some  useful  elementary  observations. 

Remarks 

1  Definiteness  of  Q  implies  spectral  stability  (i.e.,  the  spectrum  of  L  is  on 
the  imaginary  axis).  This  is  because  spectral  instability  implies  (linear  and 
nonlinear)  instability  (Lyapunov’s  Theorem),  while  definiteness  of  Q  implies 
stability  (the  Lagrange  Dirichlet  criterion). 

2  Spectral  stability  need  not  imply  stability,  even  linear  stability.  This  is 
shown  by  the  unstable  linear  system  q  =  p,p  =  0  with  a  pair  of  eigenvalues  at 
zero.  Other  resonant  examples  exhibit  similar  phenomena  with  nonzero  eigen¬ 
values. 

3  If  Q  has  odd  index  (an  odd  number  of  negative  eigenvalues),  then  L  has 
a  real  positive  eigenvalue.  This  is  a  special  case  of  theorems  of  Chetaev  [1961] 
and  Oh  [1987].  Indeed,  in  canonical  coordinates,  and  identifying  Q  with  its 
corresponding  matrix,  we  have 

L  =  JQ. 

Thus,  detL  =  det  Q  is  negative.  Since  detL  is  the  product  of  the  eigenvalues 
of  L  and  they  come  in  conjugate  pairs,  there  must  be  at  least  one  pair  of  real 
eigenvalues,  and  since  the  set  of  eigenvalues  is  invariant  under  reflection  in  the 
imaginary  axis,  there  must  be  an  odd  number  of  positive  real  eigenvalues. 

4  If  P  =  TmQ  with  the  standard  cotangent  symplectic  structure  and  if  H  is 
of  the  form  kinetic  plus  potential  so  that  an  equilibrium  has  the  form  (ge,0), 
and  if  6’2V(qe)  has  nonzero  (even  or  odd)  index,  then  again  L  must  have  real 
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eigenvalues.  This  is  because  one  can  diagonalize  62V(qe)  with  respect  to  the 
kinetic  energy  inner  product,  in  which  case  the  eigenvalues  are  evident.  In  this 
context,  note  that  there  are  no  gyroscopic  forces.  ♦ 


To  get  more  interesting  effects  than  covered  by  the  above  examples,  we 
consider  gyroscopic  systems ;  i.e.,  linear  systems  of  the  form 

Mq  +  Sq  +  Aq  =  0  (1.1) 


where  q  €  Mn,  M  is  a  positive  definite  symmetric  n  x  n  matrix,  S  is  skew,  and 
A  is  symmetric.  This  system  is  verified  to  be  Hamiltonian  with  p  =  Mq,  energy 
function 

H(q,p)=1-pM-1p+l-qAq  (1.2) 

and  the  bracket 


dF  dK 

dq‘  dpi  dq'  dpi  dpi  dpj 


(1.3) 


Systems  of  this  form  arise  from  simple  mechanical  systems  via  reduction;  this 
form  is  in  fact  the  normal  form  of  the  linearized  equations  when  one  has  an 
abelian  group.  Of  course,  one  can  also  consider  linear  systems  of  this  type  when 
gyroscopic  forces  are  added  ab  initio,  rather  than  being  derived  by  reduction. 
Such  systems  arise  in  control  theory,  for  example;  see  Bloch,  Krishnaprasad, 
Marsden,  and  Sanchez  [1991]  and  Wang  and  Krishnaprasad  [1992]. 

If  the  index  of  V  is  even  (see  Remark  3)  one  can  get  situations  where  62H  is 
indefinite  and  yet  spectrally  stable.  Roughly,  this  is  a  situation  that  is  capable 
of  undergoing  a  Hamiltonian  Hopf  bifurcation.  One  of  the  simplest  systems  in 
which  this  occurs  is  in  the  linearized  equations  about  a  special  relative  equi¬ 
librium,  called  the  “cowboy”  solution,  of  the  double  spherical  pendulum;  see 
Marsden  and  Scheurle  [1992]  and  §6  below.  Another  example  arises  from  cer¬ 
tain  solutions  of  the  heavy  top  equations  as  studied  in  Lewis,  Ratiu,  Simo  and 
Marsden  [1992].  Other  examples  are  given  in  §6.  One  of  our  first  main  results 
is  the  following: 


Theorem  1.1  Dissipation  induced  instabilities — abelian  case  Under  the 
above  conditions,  if  we  modify  (1.1)  to 

Mq  +  (S  +  tR)q  +  Ag  =  0  (1.4) 

for  small  e  >  0,  where  R  is  symmetric  and  positive  definite,  then  the  perturbed 
linearized  equations 

z  =  Lez, 

where  z  =  ( q,p )  are  spectrally  unstable,  i.e.,  at  least  one  pair  of  eigenvalues  of 
Lc  is  in  the  right  half  plane. 
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This  result  builds  on  basic  work  of  Thomson  and  Tait  [1912],  Chetaev  [1961], 
and  Hahn  [1967].  The  argument  proceeds  in  two  steps. 

Step  1  Construct  the  Chetaev  function 

W(q,p)  =  H(q,p)+0M-1p-(Aq)  (1.5) 

for  small  0  and  use  this  to  prove  Lyapunov  instability. 

This  function  has  the  key  property  that  for  0  small  enough,  W  has  the 
same  index  as  H,  yet  W  is  negative  definite ,  where  the  overdot  is  taken  in  the 
dynamics  of  (1.4).  This  is  enough  to  prove  Lyapunov  instability,  as  is  seen  by 
studying  the  equation 

W(q(T),p(T))  =  W(qo,po)  +  [T  W(q(t),p(t))dt  (1.6) 

Jo 

and  choosing  ( qo,Po )  in  the  sector  where  W  is  negative,  but  arbitrarily  close  to 
the  origin. 

Step  2  Employ  an  argument  of  Hahn  [1967]  to  show  spectral  instability. 

The  sketch  of  the  proof  of  step  2  is  as  follows.  Since  e  is  small  and  the 
original  system  is  Hamiltonian,  the  only  nontrivial  possibility  to  exclude  is  the 
case  in  which  the  unperturbed  eigenvalues  lie  on  the  imaginary  axis  at  nonzero 
values  and  that,  after  perturbation,  they  remain  on  the  imaginary  axis.  Indeed, 
they  cannot  all  move  left  by  Step  1  and  L(  cannot  have  zero  eigenvalues  since 
L€z  =  0  implies  W(z,z)  =  0.  However,  in  this  case,  Hahn  [1967]  shows  the 
existence  of  at  least  one  periodic  orbit,  which  cannot  exist  in  view  of  (1.6)  and 
the  fact  that  W  is  negative  definite.  The  details  of  these  two  steps  are  carried 
out  in  §3  and  §4. 

This  theorem  generalizes  in  two  significant  ways.  First,  it  is  valid  for  infinite 
dimensional  systems,  where  M,  5,  R  and  A  are  replaced  by  linear  operators.  One 
of  course  needs  some  technical  conditions  to  ensure  that  W  has  the  requisite 
properties  and  that  the  evolution  equations  generate  a  semi-group  on  an  appro¬ 
priate  Banach  space.  For  Step  2  one  requires,  for  example,  that  the  spectrum 
at  e  =  0  be  discrete  with  all  eigenvalues  having  finite  multiplicity.  To  apply 
this  to  nonlinear  systems  under  linearization,  one  also  needs  to  know  that  the 
nonlinear  system  satisfies  some  “principle  of  linearized  stability”;  for  example, 
it  has  a  good  invariant  manifold  theory  associated  with  it. 

The  second  generalization  is  to  systems  in  block  diagonal  form  but  with  a 
non-abelian  group.  The  system  (1.4)  is  the  form  that  block  diagonalization  gives 
with  an  abelian  symmetry  group.  For  a  non-abelian  group,  one  gets,  roughly 
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speaking,  a  system  consisting  of  (1.4)  coupled  with  a  Lie-Poisson  (generalized 
rigid  body)  system.  The  main  step  needed  in  this  case  is  a  significant  general¬ 
ization  of  the  Chetaev  function.  This  is  carried  out  in  §3. 

A  nonabelian  example  (with  the  group  SO( 3))  that  we  consider  in  §6  is  the 
rigid  body  with  internal  momentum  wheels. 

The  formulation  of  Theorem  1.1  and  its  generalizations  is  attractive  because 
of  the  interesting  conclusions  that  can  be  obtained  essentially  from  energetics 
alone.  If  one  is  willing  to  make  additional  assumptions,  then  there  is  a  formula 
giving  the  amount  by  which  simple  eigenvalues  move  off  the  imaginary  axis. 
One  version  of  this  formula,  due  to  MacKay  [1991],  states  that1 

Re\(  =  L^p^-{€  +  0(e2)  (1.7) 

where  we  write  the  linearized  equations  in  the  form 

i  =  Ltz  =  (JQ  +  cB)z.  (1.8) 

Here,  Ac  is  the  perturbed  eigenvalue  associated  with  a  simple  eigenvalue  A0  = 
iu o  on  the  imaginary  axis  at  t  =  0,  £  is  a  (complex)  eigenvector  for  Lq  with 
eigenvalue  A0,  and  (J5)anti  is  the  antisymmetric  part  of  J B. 


In  fact,  the  ratio  of  quadratic  functions  in  (1.7)  can  be  replaced  by  a  ratio 
involving  energy-like  functions  and  their  time  derivatives  including  the  energy 
itself  or  the  Chetaev  function.  To  actually  work  out  (1.7)  for  examples  like  (1.1) 
can  involve  considerable  calculation.  See  §5  for  details. 


What  follows  is  a  simple  example  in  which  one  can  carry  out  the  analysis  to 
a  large  extent  directly.  We  hasten  to  add  that  problems  like  the  double  spherical 
pendulum  are  considerably  more  complex  algebraically  and  a  direct  analysis  of 
the  eigenvalue  movement  would  not  be  so  simple. 

Consider  the  following  gyroscopic  system  (cf.  Chetaev  [1961]  and  Baillieul 
and  Levi  [1991]) 

x  —  gy  +  7X  +  ax  —  0  ) 


y  +  gx  +  6y  +  /3y  =  0 


(1.9) 


which  is  a  special  case  of  (1.4).  Assume  j  >  0  and  6  >  0.  For  7  =  6  =  0  this 
system  is  Hamiltonian  with  symplectic  form 


Q  =  dx  A  dx  +  dy  A  dy  —  gdx  A  dy  (110) 

’As  Mark  Levi  has  pointed  out  to  us,  formulae  like  (1.7)  go  back  to  Krein  [1950]  and  Krein 
also  obtained  such  formulae  for  periodic  orbits  (see  Levi  [1977],  formula  (18),  p.  33). 


7 


and  the  bracket  (1.3)  where  S  =  ^  ^  ^  and  Hamiltonian 

H  =  ^(i2+r)  +  ^(ax2+/?y2).  (Ill) 

Note  that  for  a  =  0,  angular  momentum  is  conserved  corresponding  to  the  S 1 
symmetry  of  H .  The  characteristic  polynomial  is  computed  to  be 

p(  A)  =  A4  +  (7  +  <5)  A3  +  (g2  +  a  +  0  +  y6)  A2  +  (y0  +  6a)  A  +  a/?.  (1.12) 

Let  the  characteristic  polynomial  for  the  undamped  system  be  denoted  po: 

Po(A)  =  A4  +  (<72  +  a  +  0)\2  +  a0.  (1.13) 

Since  po  is  quadratic  in  A2,  its  roots  are  easily  found.  One  gets: 

i  If  a  >  0, 0  >  0,  then  H  is  positive  definite  and  the  eigenvalues  are  on  the 
imaginary  axis;  they  are  coincident  in  a  1  :  1  resonance  for  a  =  0. 

ii  If  a  and  0  have  opposite  signs,  then  H  has  index  1  and  there  is  one 
eigenvalue  pair  on  the  real  axis  and  one  pair  on  the  imaginary  axis. 

iii  If  a  <  0  and  0  <  0  then  H  has  index  2.  Here  the  eigenvalues  may  or  may 
not  be  on  the  imaginary  axis. 

To  determine  what  happens  in  the  last  case,  let 

D  =  (g2  +  a  +  0)~  -  4 a0  =  g4  +  2 g2(a  +  0)  +  (a  -  0)2 

be  the  discriminant,  so  that  the  roots  of  (1.13)  are  given  by 

X2  =  \[-(92  +  a  +  0)±VD]. 

Thus  we  arrive  at  the  following  conclusions: 

a  If  D  <  0,  then  there  are  two  roots  in  the  right  half  plane  and  two  in  the 
left. 

b  If  D  =  0  and  g2  +  a  +  0  >  0,  there  are  coincident  roots  on  the  imaginary 

axis,  and  if  g2  +  a  +  0  <  0,  there  are  coincident  roots  on  the  real  axis. 

c  If  D  >  0  and  g2  +  a  +  0  >  0,  the  roots  are  on  the  imaginary  axis  and  if 

g2  +  a  +  0  <0,  they  are  on  the  real  axis. 

Thus  the  case  in  which  D  >  0  and  g2  +  a  +  0  >  0  (he.,  if  g2  +  a  +  0  >  2y/a0 ), 
is  one  to  which  the  dissipation  induced  instabilities  theorem  (Theorem  1.1) 
applies. 

Note  that  for  g2  +  a+0  >  0,  if  D  decreases  through  zero,  a  Hamiltonian  Hopf 
bifurcation  occurs.  For  example,  as  g  increases  and  the  eigenvalues  move  onto 
the  imaginary  axis,  one  speaks  of  the  process  as  gyroscopic  stabilization. 
Now  we  add  damping  and  get 
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Proposition  1.2  If  a  <  0,  (3  <  0,  D  >  0,  g1  +  a  +  (3  >  0  and  least  one  of  7,  8 
is  strictly  positive,  then  for  (1-9),  there  is  exactly  one  pair  of  eigenvalues  in  the 
strict  right  half  plane. 


Proof  We  use  the  Routh-Hurwitz  criterion  (see  Gantmacher  [1959,  vol.  2]), 
which  states  that  the  number  of  strict  right  half  plane  roots  of  the  polynomial 

A4  +  pi  A3  +  p2A‘  +  P3A  4-  pn 

equals  the  number  of  sign  changes  in  the  sequence 

fi  P\P2~P3  P3{P\P2  ~  Pi)  ~  PlPi  \  /in  A\ 

j  l>Pi, - > - - -  1  Pa  /  •  (114) 

l  Pi  P1P2  ~  P3  J 

For  our  case,  =  7  -f  6  >  0,  p2  —  g2  +  a  +  P  +  76  >  0, P3  =  7/?  +  a<5  <  0  and 
P4  =  aft  >  0,  so  the  sign  sequence  (1.14)  is 

{+.  +,  +.—>+}• 

Thus,  there  are  two  roots  in  the  right  half  plane.  ■ 

This  proof  confirms  the  result  of  Theorem  1.1.  It  gives  more  information, 
but  for  complex  systems,  this  method,  while  instructive,  may  be  difficult  or  im¬ 
possible  to  implement,  while  the  method  of  Theorem  1.1  is  easy  to  implement. 
One  can  also  use  methods  of  Krein  and  MacKay  to  get  the  result  of  the  above 
proposition  and  get,  in  fact,  additional  information  about  how  far  the  eigen¬ 
values  move  to  the  right  as  a  function  of  the  size  of  the  dissipation.  We  shall 
present  this  technique  in  §5.  Again,  this  technique  gives  more  specific  informa¬ 
tion,  but  is  harder  to  implement,  as  it  requires  more  hypotheses  (simplicity  of 
eigenvalues)  and  requires  one  to  compute  the  corresponding  eigenvector  of  the 
unperturbed  system,  which  may  not  be  a  simple  task. 

Example  An  instructive  special  case  of  the  system  (1.9)  is  the  system  of 
equations  describing  a  bead  in  equilibrium  at  the  center  of  a  rotating  circular 
plate  driven  with  angular  velocity  ui  and  subject  to  a  central  restoring  force. 
(These  equations  may  also  be  regarded  as  the  linearized  equations  of  motion  for 
a  rotating  spherical  pendulum  in  a  gravitational  field — see  Baillieul  and  Levi 
[1991].)  Let  x  and  y  denote  the  position  of  the  bead  in  a  rotating  coordinate 
system  fixed  in  the  plate.  The  Lagrangian  is  then 

xt  -  «y)3  +  ( yt  +  wx)2  -  ^’(z2  +  y 2)  (1-15) 

and  the  equations  of  motion  without  damping  are 

xtt  -  2uyt  +(k  -  w2)x  =  0  /,  lfi\ 

ytt  -I-  2wx(  +  (k  -  w2)x  =  0. 
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Thus  for  ui2  >  k  the  system  is  gyroscopically  stable  and  the  addition  of  Rayleigh 
damping  induces  spectral  instability.  ♦ 

It  is  interesting  to  speculate  on  the  effect  of  damping  on  the  Hamiltonian 
Hopf  bifurcation  in  view  of  these  general  results  and  in  particular,  this  example. 

For  instance,  suppose  g2  +  a  +  0  >  0  and  we  allow  D  to  increase  so  a 
Hamiltonian  Hopf  bifurcation  occurs  in  the  undamped  system.  Then  the  above 
sign  sequence  does  not  change,  so  no  bifurcation  occurs  in  the  damped  system; 
the  system  is  unstable  and  the  Hamiltonian  Hopf  bifurcation  just  enhances  this 
instability.  However,  if  we  simulate  forcing  or  control  by  allowing  one  of  7  or 
6  to  be  negative,  but  still  small,  then  the  sign  sequence  is  more  complex  and 
one  can  get,  for  example,  the  Hamiltonian  Hopf  bifurcation  breaking  up  into 
two  nearly  coincident  Hopf  bifurcations.  These  remarks  are  consistent  with  van 
Gils,  Krupa,  and  Langford  [1990]. 

The  preceding  discussion  assumes  that  the  equilibrium  of  the  original  non¬ 
linear  equation  being  linearized  is  independent  of  c.  In  general  of  course  this  is 
not  true,  but  it  can  be  dealt  with  as  follows.  Consider  the  nonlinear  equation 

x  =  f{x,e)  (1.17) 

on  a  Banach  space,  say.  Assume  /( 0, 0)  =  0  and  x(e)  is  a  curve  of  equilibria  with 
*(0)  =  0.  By  implicitly  differentiating  /(*(e),e)  =  0  we  find  that  the  linearized 
equations  at  *(c)  are  given  by 

(<5x)  =  Dx/(*(e),  e)6x 

=  Dx/(0, 0)6*  +  e  [D2f/(0,  0)6*  +  D2/(0, 0)(*'(0),  6*)] 

+0{e2)  (1.18) 

where  x'(0)  =  -Dx/(0, 0)-1/t(0, 0),  assuming  that  Dx/(0,0)  is  invertible;  i.e., 
we  are  not  at  a  bifurcation  point.  In  principle  then,  (1-17)  is  computable  in 
terms  of  data  at  (0,0)  and  our  general  theory  applies. 

A  situation  of  interest  for  KAM  theory  is  the  study  of  the  dynamics  near 
an  elliptic  fixed  point  of  a  Hamiltonian  system  with  several  degrees  of  freedom. 
The  usual  hypothesis  is  that  the  equations  linearized  about  this  fixed  point 
have  a  spectrum  that  lies  on  the  imaginary  axis  and  that  the  second  variation 
of  the  Hamiltonian  at  this  fixed  point  is  indefinite.  Our  result  says  that  these 
elliptic  fixed  points  become  spectrally  unstable  with  the  addition  of  (small) 
damping.  It  would  be  of  interest  to  investigate  the  role  of  our  result,  and 
associated  system  symmetry  breaking  results  (see,  for  example,  Guckenheimer 
and  Mahalov  [1992]),  for  these  systems  and  in  the  context  of  Hamiltonian  normal 
forms,  more  thoroughly  (see,  for  example,  Haller  [1992]).  In  particular,  the 
relation  between  the  results  here  and  the  phenomenon  of  capture  into  resonance 
would  be  of  considerable  interest. 

There  are  a  number  of  other  topics  that  should  be  investigated  in  the  fu¬ 
ture.  For  example,  the  present  results  would  be  interesting  to  apply  to  some 
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fluid  systems.  The  cases  of  interest  here,  in  which  eigenvalues  lie  on  the  imagi¬ 
nary  axis,  but  the  second  variation  of  the  relevant  energy  quantity  is  indefinite, 
occur  for  circular  rotating  liquid  drops  (Lewis,  Marsden,  and  Ratiu  [1987]  and 
Lewis  [1989]),  for  shear  flow  in  a  stratified  fluid  with  Richardson  number  be¬ 
tween  1/4  and  1  (Abarbanel  et  al.  [1986]),  in  plasma  dynamics  (Morrison  and 
Kotschenreuther  [1989],  Kandrup  [1991],  and  Kandrup  and  Morrison  [1992]), 
and  for  rotating  strings.  In  each  of  these  examples,  there  are  essential  pde  dif¬ 
ficulties  that  need  to  be  overcome,  and  we  have  written  the  present  paper  to 
adapt  to  that  situation  as  far  as  possible.  One  infinite  dimensional  example 
that  we  consider  is  the  case  of  a  rotating  rod  in  §6,  but  it  can  be  treated  by 
essentially  finite  dimensional  methods,  and  the  pde  difficulties  we  were  alluding 
to  do  not  occur.  We  also  point  out  that  some  of  the  same  effects  as  seen  here 
are  also  found  in  reversible  (but  non-Hamiltonian)  systems;  see  O’Reilly  [1993]. 
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ments  on  the  linearized  equations  in  the  block  diagonalization  theory.  We  also 
thank  the  Fields  Institute  for  providing  the  opportunity  to  meet  in  pleasant 
surroundings. 

2  The  Energy-Momentum  Method 

Our  framework  for  the  energy-momentum  method  will  be  that  of  simple  me¬ 
chanical  systems  with  symmetry.  We  choose  as  the  phase  space  P  =  TQ  or 
P  =  T’Q,  a  tangent  or  cotangent  bundle  of  a  configuration  space  Q.  Assume 
there  is  a  Riemannian  metric  ((,))  on  Q,  that  a  Lie  group  G  acts  on  Q  by 
isometries  (and  so  G  acts  symplectically  on  TQ  by  tangent  lifts  and  on  T’Q  by 
cotangent  lifts).  The  Lagrangian  is  taken  to  be  of  the  form 

L(q,v)=l-\\v\\-q-V(q),  (2.1) 

or  equivalently,  the  Hamiltonian  is 

H(q,p)=\\\p\\2q  +  V(q),  (2.2) 

where  ||  •  ||5  is  the  norm  on  TqQ  or  the  one  induced  on  T’Q,  and  where  V  is  a 
G-invariant  potential. 

With  a  slight  abuse  of  notation,  we  write  either  ( q ,  v)  or  vq  for  a  vector  based 
at  q  €  Q  and  z  =  ( q,p )  or  z  —  pq  for  a  covector  based  at  q  E  Q.  The  pairing 
between  T’Q  and  TqQ  is  written 

(pq,v  ?).  <P.v)  or  ((q,p),{q,v)}.  (2.3) 
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Other  natural  pairings  between  spaces  and  their  duals  are  also  denoted  (,). 

The  standard  momentum  map  for  simple  mechanical  G-systems  is 

3  :TQ~*  0*.  where  (J(q,v),£)  =  ((v^Q(q))) 
or  3  :  T"Q  —*  q* ,  where  (J(g,p),0  -  (p,ZQ(q))  (2.4) 

where  £q  denotes  the  infinitesimal  generator  of  f  E  0  on  Q.  We  use  the  same 
notation  for  J  regarded  as  a  map  on  either  the  cotangent  or  the  tangent  space; 
which  is  meant  will  be  clear  from  the  context.  For  future  use,  we  set  0  •  g  = 
{ZQ(q)\UQ}CTqQ. 

Assume  that  G  acts  freely  on  Q  so  we  can  regard  Q  — ►  Q/G  as  a  principal  G- 
bundle.  A  refinement  shows  that  one  really  only  needs  the  action  of  Gp  on  Q  to 
be  free  and  all  the  constructions  can  be  done  in  terms  of  the  bundle  Q  — ►  Q/G 
here,  G ^  is  the  isotropy  subgroup  for  pt  E  0*  for  the  coadjoint  action  of  G  on 
0*.  Recall  that  for  abelian  groups,  G  =  G^.  However,  we  do  the  constructions 
for  the  action  of  the  full  group  G  for  simplicity  of  exposition. 

For  each  q  E  Q,  let  the  locked  inertia  tensor  be  the  map  1(g)  :  0  — '  0* 
defined  by 

(I(?HC)  =  ((r?Q(?),<Q(?))>-  (2-5) 

Since  the  action  is  free,  1(g)  is  indeed  an  inner  product.  The  terminology  comes 
from  the  fact  that  for  coupled  rigid  or  elastic  systems,  1(g)  is  the  classical  mo¬ 
ment  of  inertia  tensor  of  the  corresponding  rigid  system.  Most  of  the  results  of 
this  paper  hold  in  the  infinite  as  well  as  the  finite  dimensional  case.  To  expe¬ 
dite  the  exposition,  we  give  many  of  the  formulae  in  coordinates  for  the  finite 
dimensional  case.  For  instance, 


Ia&  —  3%}A  aA? 

(2.6) 

where  we  write 

Kq(?)]*  =  A\(qr 

(2.7) 

relative  to  coordinates  q' , 

, i  =  1, 2, . . . ,  n  on  Q  and  a  basis  e0,  a  =  1,2,.. 

. ,  m  of 

9- 

Define  the  map  a  :  TQ  — ►  0  which  assigns  to  each  (q,v)  the  corresponding 
angular  velocity  of  the  locked  system: 

a(q,v)  =  I{q)~1(3(q,v)).  (2.8) 

In  coordinates, 

a"  =  rbgtJA\vj.  (2.9) 

The  map  (2.8)  is  a  connection  on  the  principal  G-bundle  Q  — +  Q/G.  In  other 
words,  a  is  G-equivariant  and  satisfies  a(£g(g))  =  £,  both  of  which  are  readily 
verified.  In  checking  equivariance  one  uses  invariance  of  the  metric,  equivariance 
of  J  :  TQ  — *■  0*,  and  equivariance  of  I  in  the  sense  of  a  map  I  :  Q  — *  £(0,0*) 
(j.e.,  the  space  of  linear  maps  of  0  to  0’),  namely  I (g  ■  q )  •  Adg£  =  Ad^_,I(g)  •£. 
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We  call  a  the  mechanical  connection,  as  in  Simo,  Lewis  and  Marsden 
[1991].  The  horizontal  space  of  the  connection  a  is  given  by 

hor<  =  {(?,u)  |  J(?,v)  =  0};  (2.10) 

i.e.,  the  space  orthogonal  to  the  G-orbits.  The  vertical  space  consists  of  vectors 
that  are  mapped  to  zero  under  the  projection  Q  —>  S  =  Q/G]  i.e., 

ver?  =  {£q(?)  |  £  G  0}  =  0  •  q.  (2.11) 

For  each  p  €  0*,  define  the  1-form  a ^  on  Q  by 

(a^(q),v)  =  (n,a(q,v))  (2.12) 

i.e., 

(«*)i=MW,‘-  (2-!3) 

One  sees  from  <*(£<3(5))  =  £  that  takes  values  in  J-1(p).  The  horizontal- 
vertical  decomposition  of  a  vector  ( q ,  t>)  £  TqQ  is  given  by 

v  =  hor?u  +  ver?t;  (214) 

where 

\evqv  =  [a(g,  v)]q(7)  and  hor}u  =  v  —  ver?u. 

Notice  that  hor  :  TQ  — *  J-1(0)  and  as  such,  it  may  be  regarded  as  a  velocity 
shift. 

The  amended  potential  is  defined  by 

V,(q)  =  V(q)+^(p,I(q)-^).  (2.15) 

In  coordinates, 

V,(q)  =  V(q)+±rb(q)tiapb.  (2.16) 

We  recall  from  Abraham  and  Marsden  [1978]  or  Simo,  Lewis,  and  Marsden 
[1991]  that  in  a  symplectic  manifold  (P,  Q),  a  point  ze  £  P  is  called  a  relative 
equilibrium  if 

XB(Z')eT,.(G-z.) 

i.e.,  if  the  Hamiltonian  vector  field  at  ze  points  in  the  direction  of  the  group 
orbit  through  ze.  The  Relative  Equilibrium  Theorem  states  that  if  zt  £  P  and 
ze(t)  is  the  dynamic  orbit  of  Xu  with  ze( 0)  =  ze  and  p  =  J(ze),  then  the 
following  conditions  are  equivalent 

1.  ze  is  a  relative  equilibrium 

2.  z.(t)  £  ■  ze  C  G  ■  ze 
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3.  there  is  a  £  £  g  such  that  ze(t)  =  exp(<£)  ■  ze 

4.  there  is  a  f  £  0  such  that  zt  is  a  critical  point  of  the  augmented  Hamil¬ 
tonian 

Ht(z)  :=  H(z)~  {J(z)~  n,£)  (2.17) 

5.  ze  is  a  critical  point  off/xJ  :  P  — *■  M  x  g*,  the  energy-momentum 
map 

6.  ze  is  a  critical  point  of 

7.  ze  is  a  critical  point  of  H\J~1(0),  where  O  =  G  ■  p  £  g* 


8.  [ze]  £  Pfi  is  a  critical  point  of  the  reduced  Hamiltonian  H^. 


Straightforward  algebraic  manipulation  shows  that  He  can  be 

rewritten  as 

follows 

H{(q,  v)  =  Kt-(q,v)  +  V({q)  +  (p,£) 

(2.18) 

where 

(2.19) 

and  where 

VdQ)  =  V(q)-\(^l(q)0- 

(2.20) 

These  identities  show  the  following. 

Proposition  2.1 

A  point  ze  =  ( qe,ve )  is  a  relative  equilibrium  i 

f  and  only  if 

there  is  a  £  £  g  such  that 
1-  ve  =  tQ{qe)  and 
2.  qe  is  a  critical  point  of  V( . 

The  functions  I\ %  and  are  called  the  augmented  kinetic  and  potential 
energies  respectively.  The  main  point  of  this  proposition  is  that  it  reduces  the 
job  of  finding  relative  equilibria  to  finding  critical  points  of  V 

Relative  equilibria  may  also  be  characterized  by  the  amended  potential.  One 
has  the  following  identity: 


H{q,p)  -  Kp(q,p)  +  V„(q) 

where 

AV(<7,P)  =  Tjlb-  a^(q)\\2, 
for  ( q,p )  £  3~1(p).  This  leads  to  the  following: 

Proposition  2.2  A  point  ( qe,ve )  with  3(qe,ve )  —  p  is  a  relative  equilibrium  if 
and  only  if 
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1.  ve  =  tQ{qe)  where  £  =  I  1(q )/j  and 

2.  qe  is  a  critical  point  of  V 'ii. 

Next,  we  summarize  the  energy-momentum  method  of  Simo,  Posbergh  and 
Marsden  [1990,  1991],  Simo,  Lewis  and  Marsden  [1991],  on  the  purely  La- 
grangian  side,  Lewis  [1992],  and  in  a  control  theoretic  context,  Wang  and  Kr- 
ishnaprasad  [1992].  This  is  a  technique  for  determining  the  stability  of  relative 
equilibria  and  for  putting  the  equations  of  motion  linearized  at  a  relative  equilib¬ 
rium,  into  normal  form.  This  normal  form  is  based  on  a  special  decomposition 
into  rigid  and  internal  variables. 

We  confine  ourselves  to  the  regular  case-,  that  is,  we  assume  ze  is  a  relative 
equilibrium  that  is  also  a  regular  point  (be.,  gZe  =  {0},  or  ze  has  a  discrete 
isotropy  group)  and  fx  =  3(ze)  is  a  generic  point  in  g*  (i.e.,  its  orbit  is  of 
maximal  dimension).  We  are  seeking  conditions  for  stability  of  ze  modulo  G^. 

The  energy-momentum  method  is  as  follows:  Choose  a  subspace  S  C 
kerDJ(ge ,  ve)  that  is  also  transverse  to  the  G ^  orbit  of  ( qe,ve ) 
a  find  £  £  g  such  that  8H^(ze)  =  0 
b  test  52H{(ze)  for  definiteness  on  S. 

Theorem  2.3  The  Energy-Momentum  Theorem.  If  6~H^(ze)  is  definite, 
then  ze  is  G^-orbitally  stable  in  3~1(fi)  and  G-orbitally  stable  in  P. 


For  simple  mechanical  systems,  one  way  to  choose  S  is  as  follows.  Let 


V  =  {<5g  G  TitQ  |  ((6q,XQ(qe)))  =  0  for  all  x  €  Bli¬ 
the  metric  orthogonal  complement  of  the  tangent  space  to  the  G^-orbit  in  Q. 
Let 

S  =  {6z  <E  ker  DJ(ze)  |  Tttq  ■  6z  €  V} 

where  i tq  :  T* Q  =  P  —*  Q  is  the  projection. 

If  the  energy-momentum  method  is  applied  to  mechanical  systems  with 
Hamiltonian  H  of  the  form  kinetic  energy  (A')  plus  potential  (V),  under  hy¬ 
potheses  given  below,  it  is  possible  to  choose  variables  in  a  way  that  makes  the 
determination  of  stability  conditions  sharper  and  more  computable.  In  this  set 
of  variables  (with  the  conservation  of  momentum  constraint  and  a  gauge  symme¬ 
try  constraint  imposed  on  5),  the  second  variation  of  62H ^  block  diagonalizes; 
schematically 


<52/^  = 


coadjoint  ^ 

orbit  block 

Internal  vibration 
°  block 
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Furthermore,  the  internal  vibrational  block  takes  the  form 


Internal  vibration 

'  62Vm  0 

block 

0  82I\lt  . 

where  V M  is  the  amended  potential  defined  earlier,  and  I\^  is  a  momentum 
shifted  kinetic  energy.  Thus,  formal  stability  is  equivalent  to  <52  Vj,  >0  and  that 
the  overall  structure  is  stable  when  viewed  as  a  rigid  structure,  which,  as  far 
as  stability  is  concerned,  separates  out  the  overall  rigid  body  motions  from  the 
internal  motions  of  the  system  under  consideration. 

To  define  the  rigid-internal  splitting,  we  begin  with  a  splitting  in  configura¬ 
tion  space.  Consider  (at  a  relative  equilibrium)  the  space  V  defined  above  as 
the  metric  orthogonal  complement  to  g^  •  q  in  TqQ.  Here  we  drop  the  subscript 
e  for  notational  convenience.  Then  we  split 

V  =  Vrig  ©  Vint  (2-21) 

as  follows.  Define 

vRIG  =  W?)e^liegpi}  (2.22) 

where  g^  is  the  orthogonal  complement  to  9^  in  9  with  respect  to  the  locked 
inertia  metric.  (This  choice  of  orthogonal  complement  depends  on  q,  but  we  do 
not  include  this  in  the  notation.)  From  (2.21)  it  is  clear  that  VRiG  C  V  and  that 
Vrig  has  the  dimension  of  the  coadjoint  orbit  through  //.  Next,  define 

Vint  =  {^€V|(J?,[DI(?)-^]-O=0  for  all  r 7  €  g£)  (2.23) 

where  £  =  H(g)~V.  An  equivalent  definition  is  Vint  =  {^9  G  V  |  [DI(g)_1  •  <5g]  ■ 
H  €  g^}.  The  definition  of  Vint  has  an  interesting  mechanical  interpretation  in 
terms  of  the  objectivity  of  the  centrifugal  force  in  case  G  =  50(3);  see  Simo, 
Lewis  and  Marsden  [1991]. 

Define  the  Arnold  form  A^  :  g^  x  9^  — ►  M  by 

*4/i(»7.C)  =  (a<*>iX(«,n)(C)>  =  ((*>adnX(qa »)(0)»  (2-24) 

where  \{q,n)  ■  0^  0  is  defined  by  *(?i, ,)(C)  =  I(g)-1ad^  +  ad(I(g)"V-  The 

Arnold  form  appears  in  Arnold’s  [1966]  stability  analysis  of  relative  equilibria 
in  the  special  case  Q  =  G.  At  a  relative  equilibrium,  the  form  A M  is  symmetric, 
as  is  verified  either  directly  or  by  recognizing  it  as  the  second  variation  of  Vj,  on 
Vrig  x  Vrig- 

At  a  relative  equilibrium,  the  form  A^  is  degenerate  as  a  symmetric  bilinear 
form  on  9^  when  there  is  a  non-zero  C  G  9^  such  that 

I(g)_1adJ/z  +  ad(I(g)_  V  €  9M; 

in  other  words,  when  I(9)-1  :  g*  — +  g  has  a  nontrivial  symmetry  relative  to  the 
(coadjoint,  adjoint)  action  of  g  (restricted  to  g^)  on  the  space  of  linear  maps 
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from  0*  to  g.  (When  one  is  not  at  a  relative  equilibrium,  we  say  the  Arnold  form 
is  non-degenerate  when  A^r),^)  —  0  for  ah  H  €  0^  implies  £  =  0.)  This  means, 
for  G  =  50(3)  that  A ^  is  non-degenerate  if  p  is  not  in  a  multidimensional 
eigenspace  of  I-1.  Thus,  if  the  locked  body  is  not  symmetric  (i.e.,  a  Lagrange 
top),  then  the  Arnold  form  is  non-degenerate. 

Proposition  2.4  If  the  Arnold  form  is  non-degenerate,  then 

V  =  Vr.ig  ®  Vint  •  (2.25) 

Indeed,  non-degeneracy  of  the  Arnold  form  implies  Vrig  Pi  Vint  =  {0}  and, 
at  least  in  the  finite  dimensional  case,  a  dimension  count  gives  (2.25).  In  the 
infinite  dimensional  case,  the  relevant  ellipticity  conditions  are  needed. 

The  split  (2.25)  can  now  be  used  to  induce  a  split  of  the  phase  space 

5  =  5rig©5int-  (2.26) 

Using  a  more  mechanical  viewpoint,  Simo,  Lewis  and  Marsden  [1991]  show 
how  5rig  can  be  defined  by  extending  Vrig  from  positions  to  momenta  using 
superposed  rigid  motions.  For  our  purposes,  the  important  characterization  of 
<5rig  is  via  the  mechanical  connection: 


*Srig  =  Tqatp  ■  Vrig  (2.27) 

so  5rig  is  isomorphic  to  Vrig-  Since  a M  maps  Q  to  J-1(p)  and  Vrig  C  V,  we 
get  5rig  C  S.  Define 

<Sint  =  {bz  €  S  |  6q  £  Vint};  (2.28) 

then  (2.26)  holds  if  the  Arnold  form  is  non-degenerate.  Next,  we  write 

5int  =  Wint  ©  W/NT ,  (2.29) 

where  Wint  and  W/NT  are  defined  as  follows: 

Wint  =  Tqa „  •  Vint  and  w/NT  =  {ver(7)  |  7  €  [0  ■  g]0}  (2.30) 

where  g  •  q  =  (Cq(?)  |C  €  9},  [0  ■  q]°  C  TfQ  is  its  annihilator,  and  ver(7)  £ 

TZ(T’Q)  is  the  vertical  lift  of  7  £  TfQ]  in  coordinates,  ver(g‘,7y)  =  (?' ,Pj,  0,jj). 
The  vertical  lift  is  given  intrinsically  by  taking  the  tangent  to  the  curve  <r(s)  = 
2  +  S7  at  s  =  0. 

Theorem  2.5  Block  Diagonalization  Theorem  Assume  that  the  Arnold 
form  is  nondegenerate.  Then  in  the  splittings  introduced  above  at  a  relative 
equilibrium,  62H^(ze)  and  the  symplectic  form  fi*,.  have  the  following  form: 

0  0 

62Vli  0 

0  _ 


62H((zt)  = 


Arnold 

form 

0 

0 
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and 


Q.  = 


coadjoint  orbit 

internal  rigid 

symplectic  form 

coupling 

internal  rigid 
coupling 

0 


5 

-IT 


1 

0 


where  the  columns  represent  elements  of  Srig,  Wist  and  W/NT,  respectively, 
and  I  :  W/NT  — »  Wj’NT  is  the  isomorphism  given  as  follows:  Let  wevt(y)  E  W/NT 
where  y  £  [g  ■  ^]°  and  lei  Sq  E  Vint/  then 


(I(ve rt(7)),  Tqa„  ■  Sq)  =  ( y ,  Sq) . 

As  far  as  stability  is  concerned,  we  have  the  following  consequence  of  block 
diagonalization. 

Theorem  2.6  Reduced  Energy- Momentum  Method  Let  ze  —  (qt,pe)  be 
a  ( cotangent )  relative  equilibrium  and  assume  that  the  internal  variables  are  not 
trivial;  i.e.,  Vint  ^  {0}.  If  S2H^(ze)  is  definite,  then  it  must  be  positive  definite. 
Necessary  and  sufficient  conditions  for  62 H^(ze)  to  be  positive  definite  are 

1.  the  Arnold  form  is  positive  definite  on  Vrig  and 

2.  S2Vll(qe)  is  positive  definite  on  Vint- 

This  follows  since  62Kli  is  positive  definite  and  S2H^  has  the  above  block 
diagonal  structure. 

In  examples,  it  is  this  form  of  the  energy-momentum  method  that  is  normally 
easiest  to  use. 

A  straightforward  calculation  establishes  the  useful  relation 

62VM  •  ( Sq ,  Sq)  =  S2  Vs(q')  ■  (Sq,  Sq)  +  <(D%e)  ■  Sq)f,  (Ifo.)'1  o  DI(g.)  •  Sq) 0 

(2-31) 

and  the  correction  term  is  positive.  Thus,  if  62V((qe)  is  positive  definite,  then 
so  is  S2Vft(qe),  but  not  necessarily  conversely.  Thus,  62V)i(qe)  gives  sharp  condi¬ 
tions  for  stability  (in  the  sense  of  Theorem  2.3),  while  <52Vj  gives  only  sufficient 
conditions. 

Using  the  notation  C  =  [DI-1(ge)  ■  Sq]p  E  (see  the  comments  follow¬ 
ing  (2.23)),  observe  that  the  “correcting  term”  in  (2.31)  is  given  by  (n(g« )C> C)  = 
((C<3(?e),CQ(?e)»- 

One  of  the  most  interesting  aspects  of  block  diagonalization  is  that  the  rigid- 
internal  splitting  also  brings  the  symplectic  structure  into  normal  form.  We 
already  gave  the  general  structure  of  this  and  here  we  provide  a  few  more  details. 
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We  emphasise  once  more  that  this  implies  that  the  equations  of  motion  are  also 
put  into  normal  form  and  this  is  useful  for  studying  eigenvalue  movement  for 
purposes  of  bifurcation  theory.  For  example,  for  abelian  groups,  the  linearized 
equations  of  motion  take  the  gyroscopic  form : 

Mq  +  Sq  +  Aq  =  0 

where  M  is  a  positive  definite  symmetric  matrix  (the  mass  matrix),  A  is  sym¬ 
metric  (the  potential  term)  and  S  is  skew  (the  gyroscopic,  or  magnetic  term). 
This  second  order  form  is  particularly  useful  for  finding  eigenvalues  of  the  lin¬ 
earized  equations  (see,  for  example,  Bloch,  Krishnaprasad,  Marsden  and  Ratiu 
[1991]). 

To  make  the  normal  form  of  the  symplectic  structure  explicit,  we  need  some 
preliminary  results.  See  Simo,  Lewis  and  Marsden  [1991]  for  the  proofs. 

Lemma  2.7  Let  A q  =  rjQ(qe)  £  VRIG  and  Az  =  Ta^  ■  Aq  £  5rig-  Then 

Az  =  vert  [FL(Cq(56))]  -  T'r)Q(qe)  •  pe  (2.32) 

where  £  =  I (qe)~1ad*p,  vert  denotes  the  vertical  lift,  and  F L  denotes  the  fiber 
derivative. 

Lemma  2.8  For  any  6z  £  TZeP, 

Cl{zt)(Az,6z)  =  ([DJ(ze)  •  6z],rj)  -  ((CQ(?e),  «*))•  (2-33) 

If  bz  £  Sim,  then  it  lies  in  ker  DJ,  so  we  get  the  internal-rigid  interac¬ 
tion  terms: 


Q(ze)(Az,6z)  =  -((CQ{qt),6q))  =  -  (/I.,  fo,  <*(««)]>.  (2.34) 

Since  these  involve  only  6q  and  not  bp,  there  is  a  zero  in  the  last  slot  in  the 
first  row  of  LI  and  so  we  can  define  the  operator  C  by  (2.34):  (C(Sq),  Aq)  := 
fi(zj)(Az,  8z). 

Lemma  2.9  The  rigid-rigid  terms  in  Q  are 


Lt(ze)(Aiz,A2z)  =  -(p,  [rji  ,772]),  (2.35) 

which  is  the  coadjoint  orbit  symplectic  structure. 

Next,  we  turn  to  the  magnetic  terms: 

Lemma  2.10  Let  6\z  —  Ta^  -biq  and  62z  —  Ta^-62q  £  Wint>  where  f>i q,  S2q  £ 
Vint-  Then 

Cl  (ze)(6iz,62z)  =  -dati(6iq,62q)  (2.36) 
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If  we  define  the  one  form  a ^  by  a^(q)  =  F £(£<3(9)),  then  the  definition  of 
Vint  shows  that  on  this  space  dap  =  da^ .  This  is  a  useful  remark  since  da^  is 
somewhat  easier  to  compute  in  examples.  We  also  note,  as  in  an  earlier  remark, 
that  the  magnetic  terms  can  be  equivalently  computed  from  the  magnetic  terms 
of  the  Gft  connection  rather  that  the  G  connection.  For  instance,  for  the  water 
molecule,  this  is  easier  since  in  that  case,  G  =  50(3)  while  G^  —  S1. 


Let  us  now  introduce  a  change  of  variables  r  >— *  Tttq  ■  r  of  5rig  to  Vrig  and 
p  f*  Xp  of  WjNT  to  W'[*NT,  so  that  the  representations  for  8~H^(ze)  and  QZe  will 
be  relative  to  the  space  Vrig  0  Wint  0  Wj’NT .  We  note  at  this  point  that  we 
could  have  equally  well  used  the  representation  relative  to  Vrig  0  Vint  0  Vr*NT 
and  the  results  below  would  not  materially  change  (replace  W  by  V  where 
appropriate).  Using  this  representation,  introduce  the  following  notation  for 
the  block  diagonal  form  of  82H^  : 


A»  0  0 

0  A  0 
0  0  M~l 


(2.37) 


where  A ^  is  the  co-adjoint  orbit  block;  i.e.,  the  Arnold  form,  (2  x  2  in  the  case 
of  G  =  50(3)),  A  corresponds  to  the  second  variation  of  the  amended  potential 
energy,  and  M  corresponds  to  the  metric  on  the  internal  variables. 

The  corresponding  symplectic  form  for  the  linearized  dynamics  is 


0  = 


Lm  C  0 
-CT  5  1 

0  -10 


(2.38) 


where  5  is  skew-symmetric,  1  is  the  identity  and  where  C  :  Wint  — 1 •  V^IG  is 
defined  by  (2.34).  From  the  earlier  remarks,  note  that  in  (2.37)  and  (2.38), 
the  upper  block  corresponds  to  the  “rotational”  dynamics  ( L ^  is  in  fact  the  co- 
adjoint  orbit  symplectic  form  for  G)  while  the  two  lower  blocks  correspond  to 
the  “internal”  dynamics.  In  (2.38)  C  represents  coupling  between  the  internal 
and  rotational  dynamics,  while  5  gives  the  Coriolis  or  gyroscopic  forces. 

The  corresponding  linearized  Hamiltonian  vector  field  is  then  given  by 

XN  =  (fi -1)tVF7  =  (fi -l)T62H^ 


which  a  computation  given  below  reveals  to  be 


'  -L-'Ap  0 

—L~1CM~1  ' 

r 

Xf{{r,q,p )  = 

0  0 

<1 

.  -C^L-^A^  -A 

-SM-1 

.  P  . 

(2.39) 
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where  5=5  +  CTL~1C  =  —ST. 

Thus  our  linearized  equations  have  the  form 


i  =  XH(z)  (2.40) 

where  z  =  ( r,q,p )  G  Vrig  x  Wint  x  Wt*NT  and  Xh  is  given  by  (2.39).  See  Lewis 
[1993]  and  Lewis  and  Ratiu  [1993]  for  the  explicit  expression  (in  terms  of  the 
basic  data)  of  the  linearized  equations.  To  prove  (2.39)  we  use  the  first  part  of 
the  following  lemma.  (The  remainder  of  the  lemma  will  be  used  in  our  stability 
calculations.) 


Lemma  2.11  Consider  a  vector  space  V  =  Vj  0  V2  and  a  linear  operator  M  : 
V  — *  V*  given  by  the  partitioned  matrix 


M  = 


Bn  B 12 
B21  B  22 


linear  maps.  Assume  Bn  is  an  isomorphism  and  let 

Bn  0 

3-1 


N  = 


,L  = 


1  0 

-1 


BaiBn  1  J  ’ 


d  B22 

■  v2 

— ►  V2*  are 

,P  = 

'  1 

-Bn1  Bl2 

0 

1 

0  B 22  —  B2iB11  Bi 2  J 
so  that  N  :  V  —*  V’ ,  L  :  V’  —  V*  and  P  :  V  —  V.  Then 
(i)  LMP-N 

(ii)  M  is  symmetric  if  and  only  if  Bn  and  B2 2  are  symmetric  and  Bf2  =  J52i 
(iii)  If  M  is  symmetric  so  is  N 


(iv)  If  M  is  symmetric  then  it  is  positive  definite  iff  N  is  positive  definite;  more 
generally,  the  signatures  of  M  and  N  coincide. 

Proof  (i)  is  a  computation,  while  (ii)  and  (iii)  are  obvious.  For  (iv),  note 
first  that  Lt  —  P.  If  ( ,)  denotes  the  natural  pairing,  then 

(LMPx,  x)  =  (MPx,  Ltx)  =  (MPx,Px), 

which  shows  that  N  and  M  have  the  same  signature  since  P  is  invertible.  ■ 

Lemma  2.12  The  linearized  Hamiltonian  flow  with  Hamiltonian  8~ is  given 
fey  (2.39). 

Proof  We  have  To  invert  Q,  set  Bn  =  L Bn  = 


[C  0],  B2i  = 


-CT 
0 

r- 1  - 


and  B22  = 


5  1 

-1  0 


.  From  part  (i)  of  Lemma 


2.11,  we  have  M  1  =  VN  lL  which  gives 


'  1  -L~XC  0  ' 

'  Ln  0  O' 

-1 

1  0  0' 

fi"1  = 

0  10 

0  5  1 

ctl~x  1  0 

0  0  1 

O 

1 

M 

O 

0  0  1 
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where  5  =  5  +  CT L~lC  —  —ST.  Noting  that 


we  obtain 


and 


Q-1  = 


(fi-1)7  = 


1  ■ 

-1 

'  0 

-1  ■ 

0  _ 

1 

5 

'  i;1 

0 

0 

0 

-1 

.  c7^;1 

1 

5 

_ 

0 

-L 

0 

0 

-1 

- 

;1( 

i 


Hence  we  get  result. 


3  The  Chetaev  Function  and  Lyapunov  Insta¬ 
bility 


In  this  section  we  add  a  (small)  dissipation  term  to  the  linear  Hamiltonian 
equation  (2.40)  and  show  that  this  results  in  Lyapunov  instability  for  the  linear 
system.  This  is  insufficient  to  prove  nonlinear  instability  of  the  original  system 
about  the  given  relative  equilibrium.  For  this  we  prove  a  result  on  spectral 
instability,  which  we  do  in  §4. 

We  add  dissipation  (damping)  to  the  “internal”  variables  of  the  system  only, 
in  accordance  with  the  natural  physical  models.  The  dissipation  is  assumed  to 
occur  due  to  the  addition,  to  the  Lagrangian,  of  a  Rayleigh  dissipation  function 
(see  e.g.  Whittaker  [1959]): 

K=±qTRq  =  ±(M-lp)TRM-1p,  (3.1) 

where  the  Rayleigh  dissipation  matrix  R  :  Wint  — 1 ►  Wt*NT  is  symmetric  and 
positive  definite:  R  =  RT  >  0. 

The  system  of  linearized  equations  (2.40)  becomes 


r 

4 

p 


We  note  that 


-L-'A^r-L-'CM-'p 

M~1p 

(-CTL-'A»r  -  A q-  SM~lp  -  RM~lp^j  . 


d62H{ 

dt 


-2  n. 


(3.2) 


(3.3) 
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The  presence  of  dissipation  results  in  the  addition  of  a  term  — RM~ 1  to  the 
(3,3)  block  of  the  matrix  representation  (2.39)  of  the  linear  system  (2.40). 

To  prove  Lyapunov  instability,  we  will  employ  a  generalization  of  the  Chetaev 
function  (Chetaev  [1961];  see  also  Arnold  [1987]). 


Before  doing  the  general  case,  it  is  instructive  to  analyze  the  special  case 
G  —  S1,  an  abelian  group,  where  our  system  reduces  to  the  form  of  the  system 
originally  analyzed  by  Chetaev  and  Thomson.  This  analysis  is  relevant,  for 
example,  for  examining  planar  rotating  systems  (see  e.g.  Oh  et  al.  [1989]). 

In  this  case  the  A ^  block  in  (2.37)  vanishes  and,  the  linearized  flow  X^h( 
with  the  addition  of  (internal)  damping  becomes 


q  =  M  xp 

p  =  (-A q  -  (S  +  R)M~xp) 


(3.4) 


where  R  =  RT  >  0  is  the  Rayleigh  dissipation  matrix  as  above  and  S  =  —.S'7’ 
represents  the  gyroscopic  forces  in  the  system. 

We  shall  call  (3.4)  the  Chetaev-Thomson  normal  form.  The  example 
(1.9)  analyzed  in  the  introduction  is  the  simplest  case  of  this  form. 

The  basic  question  addressed  by  Chetaev  is  the  following.  If  A  has  some 
negative  eigenvalues,  yet  the  spectrum  of 


q  =  M  xp 
p  =  -A  q-SM~lp 


is  on  the  imaginary  axis,  is  the  system  (3.4)  unstable?  Chetaev  showed  that 
this  is  indeed  the  case  for  strong  damping-,  that  is,  when  R  is  positive  definite. 
Our  proof  is  a  slight  modification  of  his.  Interestingly,  no  assumption  on  S  or 
the  the  size  of  R  is  explicitly  needed. 

Theorem  3.1  Suppose  A  is  has  one  or  more  negative  eigenvalues  and  R  is 
positive  definite.  Then  the  system  (3-4)  is  (Lyapunov)  unstable. 

The  proof  is  based  on  the  following. 

Lemma  3.2  (Lyapunov’s  Instability  Theorem)  A  linear  system  is  Lyapunov 
unstable  if  there  is  a  quadratic  function  W  whose  associated  quadratic  form  has 
at  least  one  negative  eigendirection  and  is  such  that  W  is  negative  definite. 

See,  for  example,  LaSalle  and  Lefschetz  [1963]  for  the  proof  of  this  lemma. 


To  utilize  this  lemma  to  prove  the  theorem,  we  first  assume  that  A  is  an 
isomorphism.  Let 

W(q,  p)  =  H0(q,p)  +  0Bq  ■  M~Xp,  (3.5) 
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where  Ho  is  the  Hamiltonian  for  the  undamped  system,  Hq  =  LpT M~lp  + 
\qT Aq,  0  is  a  scalar,  and  B  is  a  linear  map,  both  of  which  are  to  be  determined. 
Write 


W=-(pTqT) 


0BTM~ 1 

Calculating  the  time  derivative,  we  find  that 

W  =  -(pTqT)x 

(MT)~l  RM~l 


M-1  0(MT)~1B 

A 


(3.6) 


^((Mt)-1BA/-1  +  (Mt)-1BtM~1) 


^Bt(Mt)-1(R  +  S)M~1 


|( Mt)~1{R-S)M-1B 


^{BtM~1A  +  A(Mt)-1B) 


(3.7) 
/  P  N 
V  9  / 


Now  choose  any  positive  definite  symmetric  map  A'  :  W’int  — 1 •  Wj*NT.  Our 
choice  of  B  will  depend  on  K ,  but  K  may  be  chosen  arbitrarily,  and  this  freedom 
will  be  important  below.  We  let 


B  —  A’-1  A  :  W’int  — *■  W!*NX. 

From  Lemma  2.11,  we  see  that  W  is  negative  definite  if  and  only  if 

’  (MT)~lRM~1  +  0{(3)  0 

0  /?AA'-1  A  +  0(0-) 

is  positive  definite.  This  is  clearly  true  for  0  >  0  sufficiently  small  since 
{MT)~l RM~l  and  AK~1A  are  symmetric  and  postive  definite.  On  the  other 
hand,  by  a  similar  argument,  W  has  at  least  one  negative  eigendirection  for  0 
sufficiently  small.  Hence  by  lemma  3.2,  we  have  instability. 

To  prove  the  general  case,  in  which  A  is  allowed  to  be  degenerate,  we  proceed 
as  follows.  Split  the  space 


Wint  =  kerA  ©  (kerA)1 

into  the  direct  sum  of  the  kernel  of  A  and  its  orthogonal  complement  in  the 
inner  product  corresponding  to  M .  This  induces  a  similar  decompostion  of  the 
dual  spaces  using  M  as  an  isomorphism.  Denote  with  a  subscript  1  the  first 
component  in  this  decomposition  and  with  a  subscript  2,  the  second  component. 
In  this  decomposition,  we  have  the  block  structure 


A  = 


0  0 
0  A2  . 
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and 


M  = 


Mi  0 
0  Mo 


The  equations  (3.4)  in  this  splitting  become 


?i  =  Ml  *Pi 

qi  =  M2-1p2 

P  =  ((0,  -A252)  -  (5  +  R)M~lp). 


Notice  that  the  first  equation  for  qi  decouples  from  the  next  three  equations. 
Now  we  proceed  as  above,  with  the  function  W(q,p )  replaced  by  the  following 
function  of  (q2,p): 


W(?2 ,p)-^PTM  lp+  i?TA?2  +  /3MjK2  1A2q2  ■  (M2T)  1p2, 


where  K2  ■  (kerA)x  —  M(kerA)1)  C  Wj*NT  is  positive  definite  symmetric.  Note 
especially  that  here  we  are  using  our  freedom  to  choose  K ;  in  Chetaev,  the 
initial  choice  I\  =  A  was  made,  which  required  A  to  be  invertible.  We  now 
compute  W  as  above,  and  obtain  an  expression  similar  to  (3.7)  but  with  the 
blocks  done  according  to  the  variables  (p,q2),  and  in  which  the  top  right  and 
lower  left  expressions  are  modified,  but  are  still  multiplied  by  /?,  and  where  the 
lower  right  hand  block  is  replaced  by  the  expression  0A2K2l Ao.  Now  repeat 
the  argument  above.  ■ 


We  now  extend  our  analysis  to  the  general  equation  (2.40)  i.e.,  to  an  arbi¬ 
trary  nonabelian  symmetry  group  G.  We  show  that  indefiniteness  of  S~H^  at  a 
given  relative  equilibrium  implies  Lyapunov  instability  (again  spectral  instabil¬ 
ity  follows  from  the  analysis  in  §4). 

The  main  ingredient  is  a  generalization  of  the  Chetaev  function  (3.5).  As 
above,  we  establish  definiteness  of  the  time  derivative  of  the  function,  but  the 
analysis  is  now  more  complex.  Also  we  need  an  assumption  on  the  coupling 
matrix  C  between  the  internal  and  rotational  modes. 

Theorem  3.3  Suppose  A ^  is  nondegeneraie  and  either  A  or  .4M  has  at  least 
one  negative  eigenvalue.  Suppose  that  R  >  0  and  that  CT  is  injective.  Then  the 
system  (3.2)  is  Lyapunov  unstable. 

Proof  As  in  the  abelian  case,  we  start  with  the  assumption  that  A  is  an 
isomorphism.  In  this  case,  let 

W{q,p,r)  -  ^p-M_1p+ Ag+ ir -.4Mr 

4- 0Bq  •  M~lp  +  aDr  ■  M~lp  +  -yEr  ■  Aq  (3.8) 
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where  a,0,  and  7  are  scalars  and  B,D,  and  E  are  linear  operators,  all  to  be 
chosen.  We  write  the  matrix  representation  of  W,  in  the  ordering  (p,  q,  r)  as  : 


M~x 

w  =  \  t BtM -1 

- dtm~ 1 
L  2 


f( mt)~xb  |(. mt)~'d 

A  Ja  £ 

2 


(3.9) 


1  0  0  An  A12  A13  1  — Aj^Ait  0 

— 1  0  A22  A23  0  10 

0  0  1  J  [  Af3  Aj3  A33  J  [  0  0  1 
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0 


—  0  A22  —  A\2  —  A^A^i  A13  +  A23 

.  -^13  A\2  +  A^2  -^33 

Multiplying  (3.11)  by 

[10  -A^Al3  1 


(3.12) 


(3.13) 


on  the  right  and  by  the  transpose  of  (3.12)  on  the  left  yields 


0  A22  —  A'(2All  Au  —  Aj2A~[i  A13  +  A23 

0  —  Al{3A^  A\2  +  A23  — Af3J4uMi3  +  j433 


(3.14) 


Then  a  final  application  of  the  lemma  to  (3.13)  yields  the  block  diagonal  form 


where 


Au  0  0 

0  A22  ~  A 12-^1/ A 12  0 
0  0  2I33 


A33  =  ^33  —  A[3A^  A13  —  {Ah  —  vli3^4111v4i2)  x 

(A22  ~  A.J2A~i  Ai2)~1{A23  _  •^i2-^n1^i3)- 


(3.15) 


(3.16) 


Now  An  in  (3.10)  is  positive  definite  if  a  and  0  are  small,  since  R  is  positive 
definite.  Choose,  as  in  Theorem  3.1,  B  =  MK~lA,  and  assume  that  7  is  small, 
then  A22  —  Aj2A^Ai2  =  0AK~1A  -  Aj2A~[^ AV2.  Since  the  second  term  is  of 
higher  order  in  a,  0,y,  this  is  positive  definite.  It  remains  to  prove  positive 
definiteness  of  ^33.  Firstly,  choose 


D  =  MtK~1CtL-1A^i  :  5rig  -  Wj*N 


(3.17) 


Then  we  find: 


An  =  (MT)~l  RM~l  +  0(a)  +  0(0) 

a12  =  ?(MT)-1(R-S)K-1A-1(MT)-1CTL;1ETA 

A 13  =  j(Mt)-1(R-S)K~1CtL-1Au  +  ^M-1MK-1CtL-1A^L;1A/j 
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-  1-{MT)~XAE 

A23  =  (f  +  f )  +  ^EL-'A^ 

A22  =  0AK~xA 

A33  =  a(CTL;1AM)TK-1(CTL;1A/J). 

To  show  that  A33  is  positive  definite,  we  set  a  =  (3  =  y  and  write  it  as  a  term 
linear  in  a  plus  higher  order  terms  in  a.  Then  we  show  the  term  linear  in  a  is 
indeed  positive  definite  for  C T  injective  and  a  suitable  choice  of  E. 

We  now  isolate  the  terms  in  A33  that  are  linear  in  a.  Since 

An  =  +0(a), 

we  get  =  MR~XM(  1  +  0(a)).  Also  Ai3,A23  and  A33  are  all  O(a).  Hence 

Af3A^Ai3  =  0(a2) 

and  so  does  not  affect  definiteness,  for  small  a.  Next, 

A22  -  AJ2A^Ai2  =  aAI<~1A  +  0(a2)  =  aAA^AO  +  0(a)) 


and  so 


Also, 


(A22  -  A^A^An)'1  =  — A-1  A'A-1(1  +  0(a)). 

a 

A23  -  AJ2A±i  A13  =  A23  +  0(a“) 


and  thus 

A33  =  A33  +  0(a2)-(Aj3  +  0(a2))iA-1A'A-I(l  +  0(a))(A23  +  0(a2)) 

a 

=  A33  +  0(a2)  -  -{AT„A~X  KA~X  +  0(a2)  +  AT3A~xI< A~xO{a) 
a 

+  0(a3))(A23  +  0(a2)) 

=  A33  -  -A^A^A'A"1  A23  +  0(a2). 
a 

Hence  the  term  in  linear  in  a  is  given  by 


ch(CtL~1A^)tK~1(CtL~1A^) 


+^(aAftL-lCK-1A+^aAtlL-lETA)A-1KA-1{aAK-1CTL-1A^aAEL-lAli) 

=  (CTL-xAli)TK-x{cTL-xA^)  +  ai1l~xck~1ct  l~xa^ 

+  l-aA,L-xETCTL~lA)i  +  '-aA.L^C  EL~X  A, 

+  ^aAttL-xETKEL~lA». 
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Since  (A^L^1)  =  ~(L~1Afi)T,  the  first  two  terms  cancel  and  we  obtain 
-f  (L;'A,)t[CE  +  ( CEf  +  ^KEKL^A,). 

Now  let  E  =  -K~1Ct.  Then  CE  +  ( CE)T  +  \ET  KE  =  -fCA"1^  and 
hence 

A33  =  ^-(L;1^)T(CA-1CT)(I;1^)  +  0(a2) 

which  is  positive  definite  since  CT  is  injective  and  a  >  0. 

Since  W  is  clearly  indefinite  and  W  is  negative  definite,  we  have  Lyapunov 
instability  by  Lyapunov’s  instability  theorem. 

To  prove  the  theorem  in  the  case  that  A  is  degenerate,  split  the  variable  q 
into  ( qi ,  52)  ^  in  the  proof  of  the  abelian  case  and  note  that  the  equations  (3.2) 
decouple  into  equations  for  ?i  and  (r,  p,  52)-  Now  repeat  the  argument  using  the 
same  modifications  as  in  the  abelian  case. 

For  completeness,  we  give  the  details  in  the  extreme  case  A  =  0.  In  this  case, 
the  linearized  dynamics  with  added  dissipation  in  the  block-diagonal  normal 
form  takes  the  following  “triangular”  form: 


Mq  =  p 

p  =  -(5  +  A)A/-1p-CTA;1AMr  > 
r  =  —L~lA~1r  -  L~1CM~1p,  , 


(3.18) 


where  5  =  5  +  CT L~lC,  and  R  =  RT  >  0  is  a  matrix  of  damping  coefficients. 
Note  that  projecting  out  the  shape  variable  q  leaves  the  reduced  system  from 
(3.18)  involving  p  and  f  only,  which  can  be  handled  separately.  Let 


W(p,r)  —  ^p  M  lp  +  •  A^r  +  aDr  ■  M  lp 

where  a  is  a  scalar  and  D  is  to  be  chosen.  We  will  show  that  a  and  D  can  be  so 
chosen  that  W(p,r )  is  a  Chetaev  function  for  the  reduced  system — the  second 
and  third  equations  of  (3.18)  i.e.,  W(p,r)  is  indefinite  and  its  total  derivative 
W  along  trajectories  of  (3.18)  is  negative  definite.  This  would  then  establish  the 
Lyapunov  instability  of  the  reduced  system  and  consequently  of  the  full  system 
(3.18).  As  above,  choose 

D  =  MK-lCTL~1Ati. 

It  is  then  easy  to  verify  that 
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where, 

0n 


=  (Mt)~1RM~1  +  0(a) 


Gw  =  |(( MT)-l{R-S)K-lCTL-lA,i  +  K-1CTL-lA»L-lA „) 

022  =  a(CT  L~l  A  p)7  K~l(CT  L~l  Ap). 

By  hypothesis  Q22  >  0.  As  above,  there  is  a  range  of  a  for  which  the  matrix 


/  011  012  \ 

”  I  012  022  / 

is  positive  definite.  Further,  since  the  signature  of  a  hyperbolic  matrix  is  invari¬ 
ant  under  small  perturbations,  one  can  further  choose  a  in  the  range  (0,c)  such 
that, 


signature 


M 


-1 


a 


The  matrix 


V  2 

M-1  0 


T  A/f- 1 


DLM 


—  signature  ( 


V  0 


is  indefinite  by  hypothesis.  Thus  we  have  a  range  of 


a  for  which  W  is  a  Chetaev  function  and  we  have  proved  Lyapunov  instability. 


Remark  We  leave  it  to  the  reader  to  verify  that  standard  eigenvalue  inequal¬ 
ities  lead  to  the  condition, 

0  <  a  <  c, 


where 


c  =  min{ci,c21}, 

if  >  0 

if  Xmin(CT L~lA^lL~1C)  <  0, 

c2  =  \\{M-lRM-l)-ll2CT  L-lA^L-lC{M-lRM-l)-ll2\\ 

4- A 

max  (0f2Mf?-1M012)/A 

min  (022) 

and  ||  •  ||  denotes  the  Euclidean  norm.  ♦ 


Cl  = 


Amin  (M-'RM 


-1 


|Ami„  (CTL^ApL 


-1 


An  illustration  of  the  instability  result  of  Theorem  3.3  in  the  case  of  A  =  0, 
and  of  the  effective  use  of  the  block  diagonal  normal  form  will  be  given  in  the 
examples  in  §6. 
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4  Instability  of  Relative  Equilibria 

Our  main  result  shows  that  if  62H^  is  indefinite  at  a  given  relative  equilibrium, 
the  system  is  dissipation  unstable  about  that  equilibrium.  To  do  this,  it 
is  sufficient  to  prove  spectral  instability  of  the  linear  system  (3.2).  In  §3  we 
proved  Lyapunov  instability  of  this  system.  As  discussed  in  the  introduction, 
this  is  not  sufficient  to  prove  instability  of  the  nonlinear  system.  Hence  we 
need  to  show  that  we  do  in  fact  have  spectral  instability.  This  will  follow 
from  the  following  proposition  which  utilizes  the  eigenstructure  of  the  linearized 
Hamiltonian  system  (i.e.,  with  R  =  0)  and  a  key  observation  of  Hahn  [1967]. 

Proposition  4.1  Let  x  =  Xh(x)  be  a  linear  Hamiltonian  system.  Suppose  that 
one  adds  a  small  linear  perturbation  to  Xh  (in  particular,  a  damping  term )  and 
that  for  the  augmented  system  there  exists  a  quadratic  form  W  which  has  at  least 
one  negative  eigendirection  and  which  satisfies  W  <  0  along  the  flow.  Then  the 
augmented  system  is  spectrally  unstable. 

Proof  The  properties  of  W  imply  that  the  augmented  system  is  Lyapunov 
unstable,  as  we  have  seen.  We  now  show  that  it  is  spectrally  unstable.  Hence¬ 
forth,  we  shall  refer  to  the  augmented  system  as  the  damped  system  and  the 
perturbation  as  damping. 

We  consider  firstly  the  eigenvalue  configurations  of  the  undamped  linear 
Hamiltonian  system.  From  the  general  properties  of  Hamiltonian  matrices  (see 
e.g.  Abraham  and  Marsden  [1978])  the  possible  configurations  can  be  grouped 
into  the  following  four  categories: 

1.  There  is  at  least  one  quadruplet,  i.e.,  an  eigenvalue  configuration  shown 
in  Figure  4.1: 


Figure  4.1:  The  case  of  an  eigenvalue  quadruplet. 
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Figure  4.2:  The  case  of  eigenvalues  on  the  real  axis. 


2.  There  is  at  least  one  pair  of  real  eigenvalues,  as  in  Figure  4.2. 

3.  Neither  1.  nor  2.  holds  but  all  the  eigenvalues  are  on  the  imaginary  axis 
and  are  simple. 

4.  All  the  eigenvalues  are  on  the  imaginary  axis  and  there  is  at  least  one 
multiple  eigenvalue. 

Now  add  the  damping  terms.  In  cases  1.  and  2.,  small  damping  leaves 
eigenvalues  in  the  right  half  plane.  Hence  we  have  spectral  instability.  Now 
consider  case  3.  All  eigenvalues  cannot  move  to  the  left  half  plane  since  this 
implies  Lyapunov  stability  and  we  have  instability.  They  cannot  all  remain  on 
the  imaginary  axis  since  (for  small  damping)  they  remain  distinct  and  hence  all 
solutions  would  be  periodic  and  hence  stable.  Similarly  if  some  move  into  the 
left  half  plane  and  some  remain  (distinct)  on  the  imaginary  axis,  the  system  is 
still  stable.  The  only  remaining  possibility  is  at  least  one  moves  to  the  right 
half  plane  and  we  thus  have  spectral  instability. 

Finally  consider  case  4.  If  any  eigenvalues  move  into  the  right  half  plane  we 
have  spectral  instability  and  are  done.  Now  if  all  eigenvalues  moved  to  the  left 
half  plane  the  system  would  be  stable  and  we  know  it  is  unstable.  Similarly 
it  is  impossible  for  some  to  move  to  the  left  and  for  those  that  remain  on  the 
imaginary  axis  to  be  simple,  for  this  again  implies  stability. 

The  only  remaining  possibilities  are  a  multiple  zero  eigenvalue  or  a  multiple 
pair  of  conjugate  purely  imaginary  eigenvalues  remaining  on  the  imaginary  axis 
after  the  addition  of  damping.  We  can  show  that  both  situations  are  impossible 
for  they  contradict  W  <  0  : 

Suppose  firstly  that  there  is  a  zero  eigenvalue.  Let  W  =  zTQz  and  Xa(z)  = 
Az.  Then  W  =  zT(ATQ  +  QA)z.  But  there  exists  an  z  ±  0  and  that  Az  =  0 
and  here  W(z)  =  0,  contradicting  W  <  0. 


Now  suppose  there  is  a  pair  of  conjugate  purely  imaginary  multiple  eigenval¬ 
ues.  Then  there  exists  an  invariant  subspace  for  the  flow,  which  is  a  subspace 
of  the  generalized  eigenspace  corresponding  to  the  multiple  eigenvalues,  which 
is  invariant  for  the  matrix 

'0-60  0 

6  0  0  0 

10  0-6 
0  16  0 

Now  we  can  use  the  following  argument  of  Hahn  [1967].  There  exists  a 
solution  of  the  system  corresponding  to  (4.1)  of  the  form 

z\  —  z->  —  0  23  =  cos6t  2  4  =  sin  bt. 

However,  IT  is  a  periodic  function  of  t  when  evaluated  on  the  above  solution. 
On  the  other  hand, 

W[t)  =  W0  +  f  W(s)ds.  (4.2) 

■J  to 

Since  W  <  0,  |W|  is  bounded  away  from  zero  on  this  curve.  Hence  the  integral 
would  not  be  bounded  as  t  — ■  oo  and  W  cannot  be  periodic  on  this  trajectory. 
Hence  we  cannot  have  a  pair  of  conjugate  purely  imaginary  multiple  eigenvalues 
since  this  contradicts  W  <  0. 

Thus  we  see  that  at  least  one  eigenvalue  must  be  in  the  right  half  plane  and 
the  system  is  spectrally  unstable.  ■ 

Combining  our  results  and  using  the  notation  of  §2,  we  get 

Theorem  4.2  Assume  (for  non-abehan  groups)  CT  is  injective,  and  the  second 
variation  of  the  Energy-Momentum  function  H $  of  the  Hamiltonian  system  is 
indefinite  at  a  given  relative  equilibrium.  Then  the  addition  of  strong  ( internal ) 
Rayleigh  dissipation  gives  spectral  instability  of  the  system  about  that  relative 
equilibrium. 

The  arguments  we  have  given  are  designed  especially  to  be  applicable  to 
infinite  dimensional  systems,  even  though  we  have  so  far  confined  our  attention 
to  finite  dimensional  ones.  There  need  to  be  appropriate  assumptions  on  the 
semigroups  involved,  and  assumptions  on  the  spectra,  but  it  seems  that  the 
main  assumption  needed  for  the  above  analysis  to  be  valid  is  that  the  spectrum 
of  the  unperturbed  problem  be  discrete,  with  eigenvalues  having  at  most  finite 
multiplicity. 

Two  interesting  problems  are  the  whirling  string  and  the  rotating  circular 
liquid  drop.  We  hope  to  pursue  the  analysis  of  these  problems  using  the  present 
techniques  in  another  publication.  We  analyze  a  simple  rotating  beam,  where 
the  infinite-dimensional  calculation  reduces  to  a  finite-dimensional  one,  in  §6. 


(4.1) 
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We  now  make  some  remarks  on  the  condition  requiring  CT  to  be  injective. 


Remarks 

1  Since  CT  :  Vrig  — *•  VVfNT ,  note  it  can  be  injective  only  if  dimVRiG  = 
dimg^  <  dim Q  —  dimG,  i.e.,  dim <  d'\m(Q/G).  For  example,  for  50(3), 
this  says  that  5  <  dim  <3-  For  a  rigid  body  with  rotors  and  G  =  50(3),  this 
says  that  there  must  be  at  least  two  rotors, 

2  We  claim  that  CT  is  injective,  i.e.,  C  is  surjective,  if 

Vint  n  0  '  =  {0}- 


Proof  From  (2.34), 

(CT(Aq),8q)  =  -{(CQ(qe),6q)). 

Suppose  this  is  zero  for  all  8q  €  VinT'  Then  Ccj(?e)  €  V/nt  and  so  by  hypothesis, 
C<?(?e)  =  0-  By  freeness,  £  =  0,  and  so  by  Lemma  2.7,  £  =  I(ge)_1ad*p,  where 
T]  €  0Me,  and  so  as  rj  £  by  (2.22),  tj  =  0,  and  so  by  (2.32),  Aq  =  0.  ■ 

Notice  that  the  above  condition  is  a  hypothesis  on  Vint  being  “genuinely 
different”  from  the  “naive”  choice  of  internal  space,  namely  [g  ■  ge]1-- 

Remark  on  Internal  Symmetries:  In  some  situations,  we  will  have  internal 
symmetries  in  the  system.  In  the  Hamiltonian  case  each  such  symmetry  would 
enable  one  to  reduce  the  system  by  one  degree  of  freedom.  In  the  presence 
of  damping  (dissipation)  we  cannot  of  course  do  this,  but  one  can  nonethe¬ 
less  eliminate  the  corresponding  configuration  variables.  This  ensures  that  the 
matrix  representation  of  W  will  be  negative  semi-definite  rather  than  definite 
(with  zero  eigenvalues  due  to  the  symmetry).  The  same  analysis  as  before  then 
applies.  This  situation  will  be  illustrated  in  §6. 


5  Dissipation-Induced  Movement  of  Eigenval¬ 
ues 

In  contrast  with  the  method  of  Routh-Hurwitz  that  requires  explicit  calculations 
with  characteristic  polynomials,  the  methods  of  the  present  paper  allow  one  to 
predict  dissipation-induced  instability  of  relative  equilibria  solely  on  the  basis 
of  signature  computations — indefiniteness  of  62H^  and  injectivity  of  CT .  In  this 
sense,  the  present  paper  is  closer  in  spirit  to  the  work  of  Hermite  on  Hankel 
quadratic  forms,  cf.  the  last  chapter  of  vol  2.  of  Gantmacher  [1959].  However, 
the  classical  work  of  Routh,  Hermite  and  Hurwitz  was  aimed  at  getting  more 
refined  information — such  as  the  number  of  right  half  plane  eigenvalues — than 
just  predicting  instability.  In  the  present  context,  a  closely  related  question  is 
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that  of  determining  speeds  of  crossing  (into  the  right  half  plane)  of  pure  imagi¬ 
nary  eigenvalues  due  to  dissipative  perturbations  of  an  underlying  Hamiltonian 
system.  In  this  section  we  discuss  some  formulae  to  compute  such  speeds  and 
thereby  track  in  detail  the  mechanism  of  instability.  Our  formulae  generalize 
the  previous  work  of  Krein  [1950]  and  McKay  [1991],  and  when  specialized  to 
the  block-diagonal  normal  form  (abelian  as  well  as  non-abelian  cases)  yield  new 
and  explicit  formulae  for  crossing  speeds. 

Keeping  in  mind  the  well-known  connections  between  the  asymptotic  sta¬ 
bility  of  a  linear  system  and  solutions  to  the  matrix  Lyapunov  equation  cf. 
Bellman  [1963],  Brockett  [1970],  Taussky  [1961],  our  proofs  will  have  a  definite 
Lyapunov  theory  flavor.  In  particular,  we  will  not  need  the  Kato  perturbation 
lemma  cf.  McKay  [1991].  We  first  prove  a  basic  result. 


Lemma  5.1  Consider  a  linear  system  x  —  Ax  and  a  quadratic  form  V(x)  = 
\xTQx.  Let  V(x)  denote  the  total  derivative  ofV  along  trajectories  of  the  linear 
system,  evaluated  at  x.  Let  A  =  Ar  +  f  A,  £  spectrum(.4).  Let  f  =  xT  +  ix,  denote 
an  eigenvector  of  A  corresponding  to  A.  Then, 


Proof 


V(xr)  +  V(Xi) 
2(V(xr)  +  V(Xi))- 


(5.1) 


P(x)  =  \{*TQX  +  xTQi)  =  ]-xt(AtQ  +  QA)x,  and  thus 
V(S)  =  \?(ATQ  +  QA)t-=l-2\?Qt;  =  2\V(0. 

Thus 

.  _  v(0 
vmy 

Now 


V(*r) 

V(xr) 


Similarly, 


xJ(AtQ  +  QA)xr 
xj  Qx  r 

(XrxJ  -  A ixf)Qxr  +  xTQ( Arxr  -  A,Xi) 


xjQxr 


2Ar  -  2A, 


xjQxr 
xjQxr  ' 


V(xi) 

V(Xi) 


=  2Ar  -+•  2A,- 


xjQxr 

xjQxi 


(5.2) 


(5.3) 
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Adding  suitable  multiples  of  (5.2)  and  (5.3)  we  get, 


2Xr(xjQxr  +  xjQn) 


Hxr) 


xj  Qxr  + 


Vizi), 


V(Xr)  ^  V(Xi) 
=  2(V(xr)  +  V(xi)). 


Qxi 


Therefore, 

V(xr)+V(xi)  K(xr)  + V'(x,) 
xjQxr+xjQxi  2(V(xr)  +  V(Xi))' 


Corollary  5.2  Consider  the  matrix  Lyapunov  equation 

ATX  +  XA  =  -P  (5.4) 

associated  to  the  linear  system  x  =  Ax,  where  P  =  PT  >  0  is  given.  Suppose 
Q  =  Qt  is  a  solution  to  (5.4).  Then, 

card{A  |  A  £  spectrum(A),  Re(A)  =  Ar  >  0}  <  index(Q),  (5-5) 

where  index(Q)  means  the  number  of  negative  eigenvalues  of  Q. 

Proof  In  Lemma  5.1,  choose  Q  to  be  a  solution  to  the  Lyapunov  equation 
(5.4).  Then, 

V(x)  =  ^{xTQx  +  xTQx) 

=  l-xT{ATQ  +  QA)x 

-  —~xTPx. 

2 

From  Lemma  5.1,  for  any  eigenvalue  A  of  A, 

V(xr)  +  V(r,-) 
r  2  (V(xr)+V(Xi) 

-  ~%ixrPxr  +  XjPXj) 

( xjQxr  +  xjQxi) 

If  Pt 
2  FQS’ 

Since  P  =  PT  >  0,  Ar  >  0,  this  implies  ^ Q£  <0.  ■ 

Remark  It  is  well-known  that  when  spectrum(.4)  lies  in  the  strict  left  half 
plane,  (5.4)  has  the  unique  positive  definite  solution 
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tAT°QeAodc!. 

If  in  Corollary  5.2,  we  impose  the  additional  condition  that,  for  any  \,[i  € 
spectrum(.4),  A  +  /i  ^  0,  then,  the  inequality  (5.5)  becomes  an  equality.  This  is 
a  theorem  of  Taussky  [1961].  ♦ 

Suppose  the  linear  system  of  interest  is 


x  —  [fl  *]T Qx  +  eBx  (5-6) 

where  Q  =  — fiT  is  a  nonsingular  matrix  ( e.g .  the  symplectic  structure)  of  size 
2n  x  2n,  B  determines  a,  possibly  dissipative,  perturbation,  e  >  0  is  a  small 
parameter,  and  Q  =  QT  determines  the  energy  quadratic  form 

E(x)  =  ±xtQx  (5.7) 

for  the  underlying  unperturbed  system.  Along  trajectories  of  (5.6) 

E{x)  =  exTQBx.  (5.8) 


Corollary  5.3  Suppose  A  is  a  simple  eigenvalue  of  A  =  Q~TQ  with  eigenvector 
£  =  xr  +  ix{.  Let  A'  denote  the  real  part  of  the  eigenvalue  branch  A£  of  A(  = 
Ll~T Q  eB  emanating  from  A.  Then 


ut  U=o 

1  E'{xr)  +  E'{xi) 

2  E(xr)  +  E(xt) 


where 

E\x)  =  xtQBx. 


(5.9) 

(5.10) 


Proof  Substitute  E(x)  for  V(x)  in  Lemma  5.1  and  observe  that  simplicity 
of  A  ensures  smoothness  of  E(xer),  E(xf)  etc.  with  respect  to  e  at  e  =  0.  ■ 


If  in  Corollary  5.3,  the  eigenvalue  branch  Af  is  emanating  from  a  pure  imag¬ 
inary  eigenvalue  A  =  iw ,  then  the  formula  (5.9)  becomes  a  formula  for  the 
crossing  speed.  It  is  our  aim  to  make  this  formula  explicit  for  systems  in  block- 
diagonal  normal  form.  As  a  first  step  we  note 

Lemma  5.4  Under  the  hypotheses  of  Corollary  5.3,  and  if  A  =  iw  where  u  €  R, 


,  =  r(fl£)anti£ 


(5.11) 


where  (  •  )anti  denotes  the  anti-symmetric  part. 
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Proof 


and 

Then, 


Since  A(xr  +  ix,)  =  iu(xr  +  tx,),  we  have 
AxT  =  -UXi  =  y(£  -£), 

Axi  =  uxT  =  +0- 


E(xr)  = 


Similiarly, 


Further, 


E(x{)  = 


E'{xr)  + E'(xi) 


^xjQxr  -  ^xJnTAxr 

-yw(TQ£. 

4 


\ xjQxi  =  ixfnTJ4xi 

K^)Tn,iK+0 


xjQBxr  +  xjQBxi 
xjATnBxrAxjATQBxi 


e-e 


=  — iw£ 


T  /n5-(nfl)5 

V  2 


=  -iw£T(Q£)anti£. 


From  (5.12),  (5.13)  and  (5.14)  we  get, 

a:  = 


E'(xr)  +E'(xj)  _  |r(nBW 


(5.12) 


(5.13) 


(5.14) 


2  {E(xr)  +  E(xi))  ' 

Remark  The  special  case  fi  =  J  =  ^  o)°^  ^ormu*a  (5-H)  appears 

in  R.  McKay  [1991]  who  also  gives  it  an  averaging  interpretation.  We  note 
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that  our  result  is  a  corollary  of  the  more  general  formula  (5.9)  which  applies  to 
eigenvalues  that  are  not  necessarily  on  the  imaginary  axis.  The  proof  presented 
here  does  not  use  the  Kato  perturbation  lemma  involving  both  right  and  left 
eigenvectors — the  key  tool  in  McKay’s  argument.  ♦ 


Next  we  compute  the  average  ( E'(xr ))  over  a  cycle  of  period  2<r /ui  of  the 
periodic  solution  £e,wt  for  the  unperturbed  system.  Recall  that  at  t  =  0,  the 
formula 

E'{x  r)  =  ^(£-0TfiS(£  +  0 

holds.  For  any  other  t , 


E'{xr{t)) 


y(£e,u"  -  £e-iwt)TQB(teiut  + 1 e~iut ) 

—  {ZTQBei2ujt  -  fTQ Bi  +  £TnB(i-SrnBZe-2iult  }(5.15) 


Substituting  from  (5.15)  into  the  average  defined  by 

1  r2*/u> 

(E'(x r))  :=  —  /  E'(xr(t))dt  (5.16) 

zn/ui  J o 


we  get, 

(E'(xr))  =  l-^(fQBt,  -  FQBO  =  (5.17) 

In  evaluating  (5.16)  we  used  the  fact  that  e,ku,tdt  =  0  for  any  nonzero 
integer  k.  From  (5.14)  and  (5.11),  and  (5.12),  (5.13)  we  get, 


y  (E'(xr)) 
r  2 E(xr)  ' 


(5.18) 


This  is  the  averaging  interpretation  of  the  crossing  speed  given  by  McKay  in 
the  case  fi  =  J. 


In  the  remainder  of  this  section  we  show  how  to  adapt  the  crossing-speed 
result  (5.11)  to  the  block-diagonal  normal  form.  Recall  that  the  symplectic 
structure  of  the  block-diagonal  normal  form  is  not  canonical.  It  is  of  the  form 
“coadjoint  orbit,  internal  symplectic,  magnetic  and  coupling  terms”; 


Q  = 


C  0 
-CT  S  1 
0  -10 


The  second  variation  62H j  takes  the  form, 


Q  = 


A„  0  0 

0  A  0 
0  0  A/"1 


(5.19) 


(5.20) 
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and  the  dissipatively  perturbed  linear  system  of  interest  is  (cf.  equation  (3.2)) 
i  =  {A+eB)x  =  (n~TQ  +  eB)x  (5.21) 

where 

'00  0 

5=  0  0  0  (5.22) 

_  0  0  -RM-1  _ 

Here  x  =  ( r,q,p ),  as  in  §3.  A  key  stumbling  block  in  using  the  crossing  speed 
formula  (5.11)  is  the  need  to  calculate  the  eigenvector  £  corresponding  to  the 
eigenvalue  iw.  The  following  result  eases  the  way  a  little. 

Lemma  5.5  Consider  the  quadratic  pencil 

X2M  +  XS  +  A  -X  CT  ' 

G{  X)  =  .  (5.23) 

-XL~tC  Xl-L-JC  _ 

Then  Ao  is  a  singular  point  of  the  pencil,  i.e.,  det  [G(Ao)]  =  0,  with  corresponding 
null-vector  y0  —  ( <7o',t'o')t  iff  an  eigenvalue  of  A  with  eigenvector  £  = 

(r%,qo,Xo(Mq0)T)T. 

Proof  Note  the  equivalence  between  the  unperturbed  system  x  =  Ax  and  the 
coupled  system  consisting  of  the  second-order  internal  dynamics  together  with 
the  first  order  coadjoint  orbit  dynamics,  in  normal  form: 

M'q  +  Sq  +  Aq  =  CTr 

r  =  LfA,r+LfCq.  (5.24) 

Interpret  G( A)  as  the  Laplace  transform  representation  of  (5.24).  This  immedi¬ 
ately  identifies  singular  points  of  the  pencil  G( A)  with  the  spectrum  of  A.  The 
eigenvector-null  vector  result  is  a  direct  calculation.  ■ 

Now,  suppose  A  =  iui 0  is  a  pure  imaginary  eigenvalue  of  A  (singular  point  of 
G(A)).  Let 

<*"«)( ’”)=0' 

By  Lemma  5.5,  and  verifying  that 

.  [  0  0  0  ' 

(fiBW,  =  =  0  0  -RM-1  , 

1  [  0  RM-1  0 

we  get 

£  T(fl5)anti£  =  — *w0  qo  Rqo-  (5-25) 
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(5.26) 


Again  by  Lemma  5.5  and  (5.19), 

f  Qf  =  f0  L^ro  +  2iw0  ?0  A/<?0  +  90  S90  +  f0  Cgo  -  9o  Cr0. 

Setting, 

9o  =  V  +  */?, 
r0  =  u  4-  iv 

and  substituting  in  (5.25),  (5.26)  we  get  the  following  “block-diagonal”  version 
of  the  crossing  speed  formula, 

y _ -u)0  (rjT  Rt]  +  0T  Rf3) _ 

r  2{uT L^v  +  uT C 0  -  vT Cr,  +  rf  S0  +  u0{r,T M r,  +  0T M/3)}'  (  '  ’ 

Remark  The  crossing  speed  formula  (5.27)  can  lead  to  effective  computation 
provided  one  has  some  insight  into  det(G(f^o))  and  can  compute  a  null  vector 
of  G(iui o).  This  is  still  more  manageable  than  directly  computing  eigenvectors 
of  A  due  to  the  smaller  matrices  involved. 

Remark  In  the  abelian  case ,  =  0  =  C  and  the  crossing  speed  formula 

(5.27)  specializes  to 

y  _  -O Iq(t]T Rl]  +  PT R/3) 

r  -  2{wo(vT Mr)  +  (3TM0 )  +  r fS0)  ' 

For  a  similar  formula  for  two  degree  of  freedom  systems,  see  Haller  [1992], 
Further,  if  there  is  no  gyroscopic/magnetic  term,  i.e.,  S  =  0  then  (5.28)  predicts 
that  every  pure  imaginary  eigenvalue  of  the  unperturbed  system  is  pushed  into 
the  left  half  plane  under  a  strong  dissipation  R  >  0.  Of  course,  this  says  nothing 
about  any  eigenvalues  of  the  unperturbed  system  that  may  be  in  the  right  half 
plane — such  eigenvalues  are  bound  to  be  present  if  S  —  0  and  A  is  indefinite. 

Example  As  we  already  saw  in  the  introductory  section,  for  the  two  degrees  of 
freedom  Chetaev  problem  (c/.  equation  (1-9)), 

x  —  gy  +  67  x  +  aar  =  0 

y  +  gx+eSy  +  f3y  =  0,  (5.29) 

if  a  and  /?  are  both  negative,  and  if  we  set  e  =  0,  then  for  g2  +  a  +  0  >  2 
all  eigenvalues  are  pure  imaginary  (the  unperturbed  system  is  gyroscopically 
stable).  But,  for  strong  dissipation,  7  >  0,5  >  0  and  6  >  0,  one  pair  of 
eigenvalues  crosses  into  the  right  half  plane  and  another  pair  into  the  left  half 
plane.  This  was  shown  by  a  Routh-Hurwitz  calculation  which,  being  a  counting 
device,  is  not  capable  of  telling  us  which  eigenvalue  crosses  over  to  which  half 
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plane.  Employing  the  crossing  speed  formula  (5.28)  we  are  able  to  address 
precisely  this  problem  of  tracking  eigenvalue  movement. 

Note  that  the  Chetaev  problem  is  in  the  abelian  case  with 


One  checks  that 

det(G(iw0))  =  (-wq  +  a)(-WQ  +  P)  -  w2#2  =  0,  (5.31) 

which  determines  two  distinct  pairs  of  pure  imaginary  eigenvalues  if  g2  +  a+/3  > 
2 \/aJ5.  Corresponding  to  A  =  two,  a  null-vector  for  G(iuio)  is  given  by, 


91 

92 


igu  0 

-u)q  +  a 


t)  +  i/3. 


(5.32) 


Thus  t]  =  (0,—  ui“q  +  a)T,  (3  =  (gui0,0)T.  Substituting  in  (5.28)  we  get, 


-H^(-“o  +  a)2  +  792wn 
r~  {(-»>  +  +  s2(_w2  +  q)}- 

Using  the  relation  (5.31)  we  can  simplify  further  to  obtain 

Y  -  1  {g(-^o  +  Q')  +  7(-wg+/?)} 

2  {(—  Wq  +  a)  +  (~Wq  +  0)g2} 


(5.33) 


(5.34) 


Suppose  ±two  and  ±iu>x  are  the  distinct  eigenvalues  of  the  unperturbed  system. 
Then, 

Wg  +  w2  =  g~  +  a  +  (3. 


Therefore, 


(— w2  +  q)  +  (— w2  +  0)  +  92  —  wi  +  wo  —  ~  wo  —  wi  —  wo- 


a:  = 


Thus, 

-±{6(-u%  +  a)  +  7(-wg  +  f3)} 

K-«02} 

By  hypothesis,  a  <  0,0  <  0,  <5  >  0,7  >  0.  It  follows  that 

+  Q)  +  T(-w2  +  /?)}  >  0. 


(5.35) 


Hence  the  simple  eigenvalue  ±two  moves  to  the  right  (left)  half  plane  according 
as  whether  wj  >  wq(wi  <  wq).  See  Figure  5.1.  ♦ 
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Figure  5.1:  The  weaker  get  destabilized. 


Example  The  simplest  non-abelian  case  arises  when  G  =  50(3)  and  the  shape 
space  dimension  is  1.  A  physical  example  of  this  is  that  of  a  rigid  body  with 
an  attached  pointmass  at  the  end  of  a  spring,  free  to  oscillate  along  a  linear 
guideway.  First,  note  that  we  can  do  some  basic  calculations  without  reference 
to  a  particular  equilibrium  about  which  block-diagonal  normal  form  is  used. 
Let, 


'  0 

-9  ' 

;  C  = 

‘  Cx  ' 

■  A  — 

a  n 

“12 

.  9 

0 

.  ^2  . 

J  — 

“12 

“22 

Note  that  C7  is  not  injective,  a  case  not  covered  by  Theorem  3.3.  The  magnetic 
term  S  —  —ST  =  0,  since  the  shape  space  dimension  is  1  by  hypothesis.  Let 
A  =  a  and  M  =  m  be  the  scalar  stiffness  and  mass  respectively.  The  quadratic 
pencil  of  Lemma  5.5  takes  the  form 


r  A2m  -f  oc  —AC; 


C(A)  = 


-XC-i/g  X  - 


1 


X  Ci/g 


“n 

9 


It  can  be  verified  that 


-XC2 

“22 

9 


A  + 


“12 


p(A)  =  detG(A)  =  mA4  +  A2  ja  +  j  +  a- 


where 


Ax 


alla22  ~  al2  >  q 

9 2 
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(because  A M  >  0)  and 


^2 


C i  0^22  4-  Cofl  ii  —  2CiC2ai2 


p(C*-C.) 


a  11  <*12 
a12  a22 


>  0 


(again  because  >  0). 

There  are  three  cases  to  consider: 


(a)  If  a  >  0,  then  the  second  variation  is  postive  definite  and  all  the  roots  of 

p(A),  (he.,  eigenvalues  of  the  unperturbed  Hamiltonian  system)  are  pure 
imaginary. 

(b)  If  a  =  0,  then  there  is  repeated  root  at  the  origin  and  a  pure  imaginary 

pair  ±iu>Q. 

(c)  If  a  <  0,  two  of  the  roots  of  p(A)  are  real  with  one  root  lying  in  the  right 

half  plane. 

In  case  (a),  a  dissipative  perturbation  moves  the  eigenvalues  into  the  left 
half  plane.  This  is  already  covered  by  the  general  theory,  but  can  be  recovered 
by  the  crossing-speed  formula  (5.27)  with  5  =  0,  a  calculation  left  to  the  reader. 
Case  (c)  is  the  odd-index  case  and  the  Cartan-Chetaev-Oh  lemma  demonstrates 
instability  with  or  without  added  dissipation.  In  case  (b)  our  crossing  speed 
formula  (5.27)  applies  to  the  pair  of  pure  imaginary  roots  (since  they  are  simple). 
The  details  are  again  left  to  the  reader.  ♦ 


6  Examples 

Example  1  (The  Rigid  Body  with  Internal  Rotors)  Consider  a  rigid  body 
with  two  symmetric  rotors.  It  is  assumed  that  the  rotors  are  subject  to  a 
dissipative/frictional  torque  and  no  other  forcing.  A  steady  spin  about  the 
minor  axis  of  the  locked  inertia  tensor  ellipsoid  (be.,  the  long  axis  of  the  body), 
is  a  relative  equilibrium.  Without  friction,  this  system  can  experience  gyroscopic 
stabilization  and  the  second  variation  of  the  augmented  Hamiltonian  can  be 
indefinite.  We  aim  to  show  that  this  is  an  unstable  relative  equilibrium  with 
dissipation  added. 

The  equations  of  motion  are  (see  Krishnaprasad  [1985]  and  Bloch,  Krish- 
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naprasad,  Marsden,  and  Sanchez  de  Alvarez  [1992]): 

(Ilock  Irotor)^  —  (Ilock^  +  Irotor^r)  X  D 

fir  =  —  (Ilock  —  Irotor)  H  Ilock  ^  +  Irotor^r)  X  D  —  RQr 

(6.1) 

A  =  AQ 
0r  =  Qr. 


In  this  example,  Q  =  50(3)  x  51  x  51  and  G  =  50(3).  Also  A  €  50(3)  denotes 
the  attitude/orientation  of  the  carrier  rigid  body  relative  to  an  inertial  frame, 
Q  6  IR3  is  the  body  angular  velocity  of  the  carrier,  Qr  £  M3  is  the  vector  of 
angular  velocities  of  the  rotors  in  the  body  frame  (with  third  component  set 
equal  to  zero)  and  6r  is  the  ordered  set  of  rotor  angles  in  body  frame  (again, 
with  third  component  set  equal  to  zero).  Further,  Iiock  denotes  the  moment  of 
inertia  of  the  body  and  locked  rotors  in  the  body  frame  and  Irot0r  is  the  3x3 
diagonal  matrix  of  rotor  inertias.  We  let 

Ilock  =  diag(5i,52,53), 

Irotor  =  diag(./j,  , 0),  >  (6-2) 

Ilock  Irotor  =  diag(Ai,  A2,  A3).  ) 

Assume  that  £i  >  Bn  >  £3.  Finally,  R  =  diag(i?i ,  , 0)  is  the  matrix  of  rotor 

dissipation  coefficients,  R,  >  0. 

Consider  the  relative  equilibrium  for  (6.1)  defined  by,  fie  =  (0,0,  w)7;  ft®  = 
(0,  0, 0)T  and  6r  =  9cr  an  arbitrary  constant.  This  corresponds  to  a  steady  minor 
axis  spin  of  the  rigid  body  with  the  two  rotors  non-spinning.  Linearization  of 
the  50(3)-reduction  of  (6.1)  about  this  equilibrium  yields, 

(Ilock  ~  Irotor^ft  =  (Ilock^ft  +  Irotor^ftr)  X  fte  +  (Ijockft6)  X  <5ft 

(5Qr  =  -(Ilock  -  Irotor)-1  [(Ilock^ft  +  Irotor^r)  X 

+  (Ilock^*)  X  ft]  -  R6£lr 
89  r  =  8Qr. 

It  is  easy  to  verify  that  8CI3  =  0.  This  reflects  the  choice  of  relative  equilibrium. 
Similarly  <5ftr;s  =  0.  We  will  now  apply  Theorem  3.3  in  the  case  of  A  =  0. 

Dropping  the  kinematic  equations  for  86r  we  have  the  “reduced”  linearized 
equations 
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6Qr2 

&Q\ 

.  6O2 


-Ri 


B3  —  B2 


A\ 


LJ 


D  Bi  -  B3 

——  -u2  - - - u 

A2  A2 


0 


6  fir, 
SQr2 


0 


J  2^ 

Ai 


0 


0 


B3  -  B\ 

A3 


IjJ 


0 


6Q1 


60,2 


(6.4) 


Assume  that  lo  r£  0  (nondegeneracy  of  the  relative  equilibrium).  Then  the  above 
equations  are  easily  verified  to  be  in  the  normal  form  (3.18),  upon  making  the 
identifications,  p  =  (6Qri , 6Qr:i)T ,  q  =  (<$fii,<$Q2)T,  and, 


Since  B\  >  B3  >  B3,  A ^  is  negative  definite.  Also,  M  and  R  are  positive 
definite,  and  CT  is  injective  and  thus  all  the  hypotheses  of  Theorem  3.3  are 
satisfied.  Thus  the  linearized  system  (6.3)  or  (6.4)  displays  dissipation-induced 
instability.  ■ 


Remark  In  the  body  and  rotors  example,  the  linearized  system  was  shown 
to  be  in  block-diagonal  normal  form  by  inspection.  Our  calculations  also  re¬ 
veal  that  there  is  some  freedom  in  the  choice  of  block-diagonal  parameters-for 
instance  the  scalar  ui  could  appear  in  various  ways  in  L^,C  etc. 

Remark  This  example  is  also  instructive  in  that  we  can  verify  the  instabil¬ 
ity  result  by  a  Routh-Hurwitz  computation,  as  in  Proposition  1.2.  We  sketch 
the  computation  here  and  note  that  calculations  like  this  can  sometimes  be  te¬ 
dious,  indicating  the  usefulness  of  the  general  result,  even  in  this  relatively  low 
dimensional  case. 
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A  straightforward  calculation  yields  the  following  characteristic  polynomial 
of  the  linearized  system  (6.3)  or  (6.4),  with  no  dissipation, 


p(  A)  =  A4  +  uBA 


2.2(J|  -  (B2  -  B3))(Jl  -  (Bl  -  B3)) 


AiAi 


(6.5) 


There  are  two  eigenvalues  at  the  origin,  consistent  with  the  rank  deficit  of  2  in 

L~t  L~tC  \ 

and,  under  the  additional  physical  assumption  that 


%C) 


(Ji  -  (B2  -  B3))(J}  -  (B 1  -  B3))  >  0,  (6.6) 


the  other  two  eigenvalues  are  pure  imaginary.  In  fact,  we  assume  both  factors  in 
the  preceding  equation  (6.6)  are  negative,  since  the  rotor  inertias  are  small.  Now 
consider  the  case  in  which  Ri,R2  >  0,  and  are  small.  The  full  characteristic 
polynomial  is 


A4  +  A3(/?!  +  R2) 


2 , 2  f  R\R_ 

+  w2A2  j  + 


R: 2  ,  {Jl-iBi-BsWl-iB,  -  Bo)) ' 


A\Ao 


—  U) 


{Ri(B\  -  B3)(J|  -  (B2  -  B3))  +  R2(B2  -  B3)(J11  -  (Br  -  B3))} 


+  u 


AiA? 

2  f  (B2  -  B3)(Bx  -  B3)BlB2 
1  AiA2 


} 


(6.7) 


Now,  use  the  same  notation  for  the  characteristic  polynomial  (6.7)  as  in  Propo¬ 
sition  1.2.  We  need  to  compute  the  sign  changes  in  the  sequence  (1-14).  Clearly 
Pi  and  are  positive  since  Ri  and  R2  are  positive  and  B\  >  Bo  >  B3. 

A  computation  shows  that 

W3  o  i 

PlP2  -  P3  =  (Bl  +  R2)(RiR2)  +  — — J~J2J\ 

A\A2 

-  tV  -  B* )B2J|  +  (B2  -  Bz)j\Ri }  .  (6.8) 

AiA2 

The  first  two  terms  are  small  by  assumption  and  hence  p\p2  —  p3  is  negative.  It 
then  follows  that 

Pz{piP2  —  P3)  —  PiPa 
P1P2  ~  Pz 

is  positive. 

Hence  the  Routh-Hurwitz  sign  sequence  is  {+,+,-,+,+}  and  thus  the  ad¬ 
dition  of  dissipation  has  indeed  moved  two  eigenvalues  into  the  right  half  plane, 
causing  a  linear  instability.  ♦ 


Example  2  (Double  Spherical  Pendulum)  In  Marsden  and  Scheurle  [1992] 
the  double  spherical  pendulum  is  discussed.  In  particular,  relative  equilibria, 
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called  “cowboy  solutions”  are  found  explicitly  and  have  a  shape  in  which  the 
horizontal  projections  of  the  two  rods  point  in  opposite  directions.  The  group 
in  this  case  is  S1 ,  corresponding  to  rotations  about  the  vertical  axis.  It  is 
verified  that  indeed  the  linearized  equations  are  in  our  standard  form  Mq  + 
Sq  +  \q  =  0,  but  where  the  3x3  matrices  M,  S  and  A  have  extra  zeros  due 
to  discrete  symmetries.  It  is  found  that  in  large  regions  of  parameter  space 
(determined  by  the  pendulum  lengths,  masses  and  angular  momentum),  that 
A  has  signature  while  the  eigenvalues  of  the  linearized  system  lie  on 

the  imaginary  axis.  It  follows  from  Theorem  1.1  or  4.2  that  if  one  adds  joint 
friction  (so  that  the  total  angular  momentum  is  still  conserved)  then  the  cowboy 
solutions  become  spectrally  unstable.  This  example  is  a  good  one  in  that  direct 
analytical  computation  of  eigenvalue  movement  to  see  this  instability  would  be 
quite  complicated.  We  also  point  out  that  experiments  of  John  Baillieul  (Boston 
University)  confirm  this  instability.  ■ 

We  also  point  out  that  similar  eigenvalue  and  energetic  situations  arise  in  a 
number  of  other  examples;  among  them  are: 

1.  The  heavy  top —  see  Lewis,  Ratiu,  Simo  and  Marsden  [1992] 

2.  The  rotating  liquid  drop—  see  Lewis  [1989] 

3.  Shear  flow  in  a  stratified  fluid  with  Richardson  number  between  1/4  and 
1;  see  Abarbanel  et  al.  [1986] 

4.  Plasma  dynamics;  see  Morrison  and  Kotschenreuther  [1989],  Kandrup 
[1991],  and  Kandrup  and  Morrison  [1992], 

The  last  three  examples  mentioned  are  infinite  dimensional,  which  provide  mo¬ 
tivation  for  extending  our  methods  to  cover  such  cases.  One  infinite  dimensional 
example  we  can  handle  is  the  next  one. 

Example  3  We  now  consider  a  partial  differential  equation  for  which  one  can 
analyze  dissipation  induced  instability  by  finite-dimensional  techniques.  We 
consider  a  Lagrangian  for  a  model  of  a  nonplanar  rotating  beam  with  “square” 
cross-section.  The  beam  is  assumed  to  be  of  Euler-Bernoulli  type.  It  is  fixed 
to  the  center  of  a  circular  plate  rotating  with  constant  angular  velocity  ui,  with 
undeflected  position  perpendicular  to  the  plate  along  the  z-axis  of  a  Cartesian 
coordinate  system  fixed  in  the  plate.  The  beam  is  inextensible  and  can  deflect 
in  the  x-  and  y-directions.  (The  planar  version  of  this  model  is  analyzed  in 
Baillieul  and  Levi  [1987].)  The  Lagrangian  is  chosen  to  be 

L{x,y,xt,yt)  =  -j  ({xt  -  yu)2  +  (yt  +  xu)2)  dz 

~\J  kiy'zz+x2z)d:.  (6.9) 
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where  k  is  an  elastic  constant. 

The  equations  of  motion  with  Rayleigh  damping  and  damping  constant  7 
are: 

xtt  -  2 uyt  -  w2x  +  kxzzzz  +  yxt  =  0 
ytt  +  2 uxt  -  ui2y  +  kyzzzz  +  7 yt  =  0. 

The  natural  boundary  conditions  are: 

x(0)  =  x'(0)  =  x"(l)  =  x"'(l)  =  0 

2/(0)  =  2/(0)  =  2/"(l)  =  2/"'(l)  =  0 

where  '  denotes  the  z-derivative. 

The  equilibrium  states  are  given  by  xn  —  ytt  =  0,x(  =  yt  =  0  and  hence 

-w2x  +  kxZZZ2  =  0 

-ui2y  +  kyZZZ2  =  0. 


(6.10) 


(6.11) 


(6.12) 


We  set  k  =  1  for  convenience. 

The  fourth  order  operator  with  the  given  boundary  conditions  has  compact 
inverse  and  hence  (see  e.g.  Baillieul  and  Levi  [1987])  the  eigenvalues  of  equations 
(6.9)  are  given  by  0  <  u2  <  w2  <  . . .  — ►  00  with  corresponding  eigenfunctions 
x(z)  =  x,(z)  and  y(z)  =  p,(z)  respectively.  By  our  choice  of  the  elastic  con¬ 
stants,  x,(z)  =  j/i(z). 

Consider  now  the  undamped  case,  7  =  0,  and  write  the  equations  in  first 
order  form,  letting  qx  —  x,  92  =  2/1P1  =  x,,p2  =  2/<-  We  obtain: 


Qu  =  Pi 
92 1  =  P2 

Pit  =  2uip2+w2qx  -  qUzZ2 
P2t  =  —2utpi  +ui2qn  —  qZzzzz- 

Let  z  —  [?!  g2  Pi  P2]t-  The  system  is  thus  of  the  form  z<  =  Az  where 


(6.13) 


0  0  10 
0  0  0  1 
-<94+w2  0  0  2w 

0  -34+w2  -2w  0 


(6.14) 


The  stability  of  the  equilibria  are  determined  by  the  eigenvalues  of  A.  In  ad¬ 
dition  to  the  zero  eigenvalue,  one  can  check  that  A  has  eigenvalues  ±i(u>  ±uj) 
with  corresponding  eigenvectors 


xi 

’  xi 

iXj 

—  ixj 

i(ui  +u>j)xj 

-i(ui  +U!j)Xj 

— i(u  +  uij)xj 

_-(w  +  Wi)x>  . 

'  XJ 

ixj 

—  ixj 

i(ui  ~  u)xj 

—  i(u>j  —u>)xj 

.  -(^;  -U>)Xj  _ 

m-(uj-u)xj  _ 
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Now  project  the  system  onto  the  invariant  subspace  spanned  by  the  four 
eigenvectors  corresponding  to  the  eigenvalues  ±i(w  ±  uj).  We  see  that  on  this 
subspace  we  have  a  gyroscopic  system  in  Chetaev  normal  form  (1.9).  In  fact,  it 
is  identical  to  the  system  describing  the  rotating  bead  given  in  §1,  with  spring 
constant  k  =  wj.  Hence  for  w 2  >  wf  we  can  see  that  addition  of  dissipation 
causes  the  system  to  become  spectrally  unstable.  In  fact,  for  wj  <  w2  <  wj  +  1 
there  are  j  gyroscopically  stable  Chetaev  subsystems  whose  eigenvalues  will  be 
driven  into  the  right-half-plane  on  the  addition  of  dissipation.  ■ 
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